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NUMERICAL  AND  ALGEBRAIC 


Gain  in  decibels  of  P2  relative  to  Pi 

G  =  101og10(P2/P1). 


To  within  two  percent 

(2t r)1/2  «  2.5;  7 r2  «  10;  e3  «  20;  210  «  103. 


Euler-Mascheroni  constant1  7 

=  0.57722 

Gamma  Function  T(x  +  1)  =  x 

T(x): 

r(l/6)  =  5.5663 

r(3/5)  =  1.4892 

T(l/5)  =  4.5908 

T(2/3)  =  1.3541 

r(l/4)  =  3.6256 

T(3/4)  =  1.2254 

T(l/3)  =  2.6789 

T(4/5)  =  1.1642 

T(2/5)  =  2.2182 

T(5/6)  =  1.1288 

T(l/2)  =  1.7725  =  v/7 

r(i)  =  1.0 

Binomial  Theorem  (good  for  |x|  <  1  or  a  —  positive  integer): 


(1  -f-  x)°  =  ^  ^  Qik  =  l  +  ai  + 


«(«  -  !)  2  ,  «(«  “  !)(«  -  2)  3 


2! 


x  + 


3! 


x  +  .  . 


k  =  0 


Rothe-Hagen  identity2  (good  for  all  complex  x,  y,  z  except  when  singular) 


Ex  /X  +  kz\  y  / 

x  +  kz  '  k  '  y  +  (n  —  fc)z  ' 


y  +  (n-  A:)z> 
n  —  k 


k  =  0 


z  +  y  /x  +  y  +  71Z\ 
x  +  3/  4-  nz  \  n  ' 


Newberger’s  summation  formula3  [good  for  fi  nonintegral,  Re  (ot  +  p)  >  -  1  ] 

OO 

(-1)"  Ja-yn(z)J0+yfl(z)  7 


E 


n  —  —  oc 


n  +  /j. 


Sin  (ITT 


Jo  4-  -r  (  z  )  J  3  — 1  m  (  z  )  • 
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VECTOR  IDENTITIES4 

Notation:  /,  g,  are  scalars;  A,  B,  etc.,  are  vectors;  T  is  a  tensor;  /  is  the 


unit  dyad. 

(1) 

A  BxC  = 

:  AxB-C  =  B  CxA  =  Bx< 

CA  : 

=  C 

AxB  = 

Cx A  B 

(2) 

A  x  (B  x 

C)  =  (C  x  B)  x  A  =  (A  •  C)B  - 

-(A 

•  B)C 

(3) 

X 

X 

< 

C)  +  B  x  (C  x  A)  +  C  x  (. 

A  x  B)  = 

0 

(4) 

(A  x  B)  • 

(C  x  D)  =  (A  ■  C)(B  •  D)  - 

-(A- 

D)( 

B  C) 

(5) 

(A  x  B)  x  (C  x  D)  =  (A  x  B  •  D)C  - 

-(A 

x  B 

•  C)D 

(6) 

V(/<7)  =  V(<//)  =  fVg+gVf 

(7) 

V  •  (/A)  = 

=  /V  -  A+  A  •  V/ 

(8) 

V  x  (/A) 

=  /V  x  A  +  V/  x  A 

(9) 

V  ■  (A  x  1 

B)  =  B  VxA-A  VxB 

1 

(10) 

V  x  (A  x 

B)  =  A(  V  •  B)  -  B(  V  •  A) 

+  (B 

•  V)A  -  (A  • 

V)B 

(11) 

A  x  (V  x 

B)  =  (VB)  •  A  -  (A  •  V)B 

(12) 

V( A  •  B) 

=  A  x  (V  x  B)  +  B  x  (V  x 

A)  + 

(A- 

V)B  +  ( 

;b  •  v)a 

(13) 

V2/  =  V 

-  V/ 

(14) 

V2A  =  V(V  -  A)  -  V  x  V  x  A 

(15) 

V  x  V/  = 

:  0 

(16) 

V  •  V  x  A 

.  =  0 

If  ei,  e2,  e3  are 

orthonormal  unit  vectors,  a 

second-order  tensor  T can 

be  written  in  the  dyadic  form 

(17)  T  —  .  T{  j  Gi  &j 

In  cartesian  coordinates  the  divergence  of  a  tensor  is  a  vector  with  com¬ 
ponents 

(18)  (V.T),  = 

[This  definition  is  required  for  consistency  with  Eq.  (29)].  In  general 

(19)  V  •  (AB)  =  (V  •  A)B  +  (A  ■  V)B 

(20)  V  •  (fT)  —  V/T  +  /V-T 
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Let  r  =  ix  4- }y  +  kz  be  the  radius  vector  of  magnitude  r,  from  the  origin 
to  the  point  x ,  y,  z.  Then 

(21)  V  •  r  =  3 

(22)  V  x  r  =  0 

(23)  Vr  =  r/r 

(24)  V(  1/r)  =  -r/r3 

(25)  V  •  (r/r3)  =  47r<5(r) 

(26)  Vr  =  / 

If  V  is  a  volume  enclosed  by  a  surface  S  and  dS  =  ndS,  where  n  is  the 
unit  normal  outward  from  V, 

(27)  [  dWf  =  [  dSf 


(28)  /  dVV  A  =  /  dS  ■  A 


-I, 


(29)  /  dW-T  =  dS-T 


-I, 


(30)  /  x  A  =  dS  x  A 


(31)  /  dV(fV2g  -  gV2f)  =  /  dS  ■  (fV g  -  gV f) 


(32)  /  dV  (A  VxVxB  —  B-VxVxA) 


<£S-(BxV  xA-AxVxB) 


If  S  is  an  open  surface  bounded  by  the  contour  C,  of  which  the  line 
element  is  dl, 


J  dS  x  V/  =  <j) 
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/, 


(34)  j  <iS  •  V  x  A  =  (h  d  1  •  A 

s  J  c 


(35) 


J  ( dS  x  V)  x 

i 


A  =  (b  d  1  x  A 

c 


(36)  j  dS  •  (V/  x  V,?)  =  j)  fdg=  -  J)  gdf 


DIFFERENTIAL  OPERATORS  IN 
CURVILINEAR  COORDINATES5 

Cylindrical  Coordinates 

Divergence 

__  .  1  d  f  A  x  1  aA*  3AZ 

V  ■  A  =  —  — — (r  Ar)  + - 1-  - 

r  dr  r  d<(>  dz 

Gradient 

(V/)r=  |£;  (V/),  = 

or  r  d<(>  dz 


Curl 


(V  x  A)r  = 


1 aAz  a  a 

r  d4>  dz 


<t> 


(V  x  A)*  = 


aAr  3A, 
dz  dr 


Id  1  aAr 

(VxA),  =  -— (rA^)---f 
r  dr  r  d<p 


Laplacian 


V2/ 


1  d_  (  d£\  1  d2f  d2f 

r  dr  \  dr  j  r 2  d<p 2  az2 
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Laplacian  of  a  vector 

_  2  2  dA 4,  A r 

<V  A  )r  =  v*Ar--—-- 


(V2A)*  =  V 2  A*  + 


2  dAr  A $ 

r2  d(f>  r2 


(  V2  A)Z  =  V2  A; 


Components  of  (A  •  V)B 

dBr  A *  dBr  ,  .  dBr  A+B+ 

(A  •  VB)r  =  Ar  — —  +  —  — -  +  Az— - 

v  dr  r  d(j>  dz  r 


dB *  A+  dB +  ,  4  dBj,  A4,Br 
(A.VB  U  =  Ar-±  +  -±^±  +  A,—  +  — 


dB,  A#  dB,  .  8B, 
(A.VB),  =  Ar—  +  —  _  +  A,— 


Divergence  of  a  tensor 

(x?  T\  —  -  d  /  T  x  ,  1  dT^, r  dTz r  _  r^'<j>4> 
(V  '  7)r  “  r  Trr)  +  r  d<t>  +  dz 


(v  ■ T>*  = ;  57(rr-»)  + ;  +  “ar  +  ~r 


_  1  d  ,  m  ,  1  dT*z  dTzz 
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Spherical  Coordinates 


Divergence 


1  d 


1  d 


1  dA, 


V  •  A  =  --(rMr)+  — ; - —  (sin  0,4*)  +  — - ^ 

r2  dr  r  sin  9  dO  r  sin  0  dtp 


Gradient 


<v/)r  =  y-,  (v/)s  =  (v/>*  =  — - y 

or  t  dO  t  sin  0  dtp 


Curl 


(V  x  A)r  = 


1  d  ,  ■  ,  1  »A, 

- (sin  B Aa.)  — - - — 

r  sin  0  dO  r  sin  0  d<f> 


^  1  dAr  1  9  , 

(V  x  A)e  =  — —  — - — —  (t A <f>) 

r  sin  0  dtp  r  or 


Id  1  dAr 

(VxA  U  =  -7r(rA,) - -f- 

r  dr  r  dO 


Laplacian 


V2/ 


=  LI l  (r2^  +  1  d  (s  e—}  + - 

r2  dr  \  dr  )  r2  sin  0  dO  \  dO  )  r2  sii 


_1 _ d2f 

sin2  6  dtp2 


Laplacian  of  a  vector 


( V2 A)r  =  V2Ar  -  2-L 

rz 


2  dAg  2  cot  OAff  2  dA * 

r 2  dO  r2  r 2  sin  0  dtp 


(V"A)e  =  +  — 


2  dA, 


2  cos  0  dA$ 


t 2  dO  r 2  sin2  9  r2  sin2  6  dtp 


2  2  Ad>  2  dAr  2  cos  0  dAe 

(V2A)*  =  v2a*  -  --V-  +  r  —  +  -2  -2-c^n- 

rl  smz  0  sin  0  dtp  rz  sinz  a  dtp 
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Components  of  (A  •  V)B 


(A  •  VB)r  =  Ar 


dPr 

dr 


Ae  dDr 

r  ae 


A 4  dPr 
rsin0  d<p 


Ae  P$  +  A$  P<p 

r 


(A  •  VB)„  =  Ar 


dPe 

dr 


Ae  dPe 

~  ~w 


+ 


A  4,  dP0 

r  sin  B  d<f> 


+ 


Ae  Pr 

r 


cot  BA^P^ 


T 


dPj,  Ae  dP 4  Aj,  dP 4  AtPr  cot  BA^Pg 

(A  •  VB)*  =  Ar—^  + - r—  +  — r-r—  +  -  +  - - 


dr 


t  dB  r  sin  B  d<p 


Divergence  of  a  tensor 

_  1  d 

(V  .  T)r  = 

r2  dr 


(V-  T)e  =  ^^(r2Tr-e)  + 


(V  -  T)+  =  \~{r2Tr4>)  4- 

r2  dr  r  sin  B  dB 


r  sin  8  dB^ 

+  ; 

1 

s, 

r  sin  8  dBK 

+ 

1 

d 

1  dT^r  Tee  +  T 4,4, 


rsinB  d<f> 


r  sin  B  d<f> 
(sin  BT0<t>) 


dT^e  Te, 

+  — 


cot  BT^t 


1  dT^t  T^r  cot  BT^e 
+ - —  +  -2—  +  - — - 

r  sin  B  d4>  r  r 
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DIMENSIONS  AND  UNITS 


To  get  the  value  of  a  quantity  in  Gaussian  units,  multiply  the  value  ex¬ 
pressed  in  SI  units  by  the  conversion  factor.  Multiples  of  3  in  the  conversion 
factors  result  from  approximating  the  speed  of  light  c  =  2.9979  x  1010  cm/sec 
~  3  x  1010  cm/sec. 


Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

Capacitance 

C 

t  q 

l 

farad 

9  x  1011 

cm 

ml 2 

Charge 

coulomb 

3  x  109 

statcoulomb 

1 1 

<7 

t 

Charge 

density 

Q 

m1 /2 

coulomb 

/m3 

3  x  103 

statcoulomb 
/  cm3 

P 

l3 

l3/2t 

Conductance 

tq 2 

l 

siemens 

9  x  1011 

cm/sec 

ml2 

t 

Conductivity 

a 

tq 2 

1 

siemens 

9  x  109 

sec  - 1 

ml3 

t 

/m 

Current 

/,  i 

q 

m1/2/3/2 

3  x  109 

statampere 

ampere 

t 

t2 

1/2 

Current 

density 

q 

m  ' 

t  ; 

ampere 

/m2 

3  x  105 

statampere 
/  cm 

J  *  J 

r-t 

ll/2t2 

rn 

m 

g/cm3 

Density 

p 

l3 

l3 

kg/m3 

10"3 

m1'2 

Displacement 

D 

H 

coulomb 

12? r  x  105 

statcoulomb 

i 2 

D/2t 

/m2 

/  2 
/cm 

Electric  field 

E 

Till 

m1/2 

volt/ m 

1  -4 

-  x  10 

3 

statvolt/cm 

t2  q 

ll/2t 

Electro- 

nil 2 

rnl/2il/2 

volt 

1 

-  x  10  " 

statvolt 

motance 

Em  1 

t2q 

t 

3 

ml2 

rnl 2 

Energy 

J.  IV 

t2 

t2 

joule 

10' 

erg 

Energy 

density 

W.  f 

in 

u 2 

in 

It2 

joule/m3 

10 

erg/cm 

10 


Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

hoi 

SI 

Gaussian 

II  nits 

Factor 

U  nits 

ml 

ml 

orcc 

F 

newton 

105 

dyne 

1 

1 

hertz 

requency 

t 

t 

1 

hertz 

ml 2 

t 

1  ii 

npedance 

Z 

- 

ohm 

-  x  10  11 

sec  /  cm 

tq 2 

1 

9 

iductance 

ml 2 

t 2 

1  ii 

L 

<72 

l 

henry 

-  x  10  11 

9 

sec2/ cm 

cngth 

l 

l 

l 

meter  (m) 

102 

centimeter 

(cm) 

m1/2 

lagnctic 

intensity 

H 

ampere- 

turn/m 

4tt  x  10~3 

oersted 

It 

P/2  t 

lagnetic  flux 

<I> 

ml 2 

ml'2l3''7 

108 

maxwell 

weber 

<7 

m 

t 

1/2 

TU  ' 

lagnetic 

B 

tesla 

104 

gauss 

tq 

iim 

induction 

lagnetic 

m,  ft 

l2q 

m1/2!5/2 

ampere-m2 

103 

oersted- 

t 

moment 

t 

cm3 

1/2 

lagnetization 

M 

<1 

771  7 

ampere- 
turn  /m 

10“3 

oersted 

It 

W2t 

lagneto- 

M, 

<7 

m1/2/1/2 

47T 

gilbert 

ampere- 

t 

t2 

10 

motance 

Mmf 

turn 

lass 

m,  M 

m 

771 

kilogram 

103 

gram  (g) 

(kg) 

lomentum 

ml 

ml 

kg-m/s 

g  cm / sec 

pP 

t 

t 

105 

771 

771 

10-1 

lomcntum 

Th 

kg/m2-s 

g/cm2  sec 

density 

ml 

henry/ m 

1  7 

ermeability 

— 

i 

—  x  10 

72 

_ 

47T 

11 


m 

Dimensions 

Physical 

9 

SI 

Conversion 

Gaussian 

Quantity 

SI 

Gaussian 

Units 

Factor 

Units 

Permittivity 

€ 

#2„2 
t  q 

ml3 

1 

farad/m 

36tt  x  109 

— 

1/2 

Polarization 

9 

m  ' 

coulomb/m2 

3  x  105 

statcoulomb 
/  2 

a  ; 

l2 

l'/2t 

/crn 

Potential 

ml2 

m ‘/’f1'2 

volt 

1  -2 
-  x  10  2 

statvolt 

t2q 

t 

3 

ml2 

ml2 

Power 

t3 

t3 

watt 

107 

erg/sec 

Power 

density 

m 

m 

Ti3 

li3 

watt/m3 

10 

erg/cm3-sec 

Pressure 

■ 

771 

m 

dyne/cm2 

mfiif 

li2 

It2 

pascal 

10 

Reluctance 

n 

92 

1 

ampere-turn 

4x  x 10-9 

cm-1 

Till2 

7 

/weber 

ml~ 

t 

-  x  10-11 

Resistance 

R 

ohm 

sec/cm 

tq2 

1 

9 

ml3 

1  9 

Resistivity 

V  ,P 

tq2 

t 

ohm  m 

-  x  10-9 

9 

sec 

Thermal  con- 

k,  k 

ml 

ml 

watt/m- 

10s 

erg/cin-sec- 

ductivity 

t3 

t3 

deg  (K) 

deg  (K) 

Time 

t 

t 

t 

second  (s) 

1 

second  (sec) 

ml 

Vector 

A 

weber/m 

106 

gauss-cm 

tq 

potential 

t 

l 

l 

Velocity 

V 

t 

t 

m/s 

102 

cm/sec 

m 

m 

kg/ m-s 

Viscosity 

It 

It 

10 

poise 

1 

1 

Vorticity 

C 

t 

t 

s-  1 

1 

sec-1 

ml2 

ml 2 

Work 

w 

t2 

t2 

— 

joule 

107 

erg 
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INTERNATIONAL  SYSTEM  (SI)  NOMENCLATURE6 


Physical 

Name 

Symbol 

!  Physical 

Name 

Symbol 

Quantity 

of  Unit 

for  Unit 

Quantity 

of  Unit 

for  Unit 

^length 

meter 

m 

electric 

potential 

r 

volt 

V 

*mass 

kilogram 

kg 

electric 

ohm 

Q 

*time 

second 

s 

resistance 

^current 

ampere 

A 

electric 

conductance 

siemens 

S 

’•‘temperature 

kelvin 

K 

electric 

farad 

F 

•amount  of 

mole 

mol 

capacitance 

substance 

magnetic  flux 

weber 

Wb 

^luminous 

candela 

cd 

intensity 

magnetic 

inductance 

henry 

H 

fplane  angle 

radian 

rad 

magnetic 

tesla 

T 

fsolid  angle 

steradian 

sr 

intensity 

frequency 

hertz 

Hz 

luminous  flux 

lumen 

lm 

energy 

joule 

J 

illuminance 

lux 

lx 

force 

newton 

N 

activity  (of  a 
radioactive 

becquerel 

Bq 

pressure 

pascal 

Pa 

source) 

power 

watt 

W 

absorbed  dose 
(of  ionizing 

gray 

Gy 

electric  charge 

coulomb 

C 

radiation) 

*SI  base  unit  fSupplementary  unit 


METRIC  PREFIXES 


Multiple 

Prefix 

Symbol 

Multiple 

Prefix 

Symbol 

10“l 

MEM 

10 

deca 

da 

10~2 

WBm 

102 

hecto 

h 

10~3 

— 

103 

kilo 

k 

10“6 

micro 

106 

mega 

M 

10~9 

nano 

n 

109 

giga 

G 

10-12 

pico 

1()12 

tera 

T 

10-15 

fern  to 

■ 

1015 

peta 

P 

10-18 

atto 

10IS 

exa 

E 

13 


PHYSICAL  CONSTANTS  (SI)7 


Physical  Quantity 

Symbol 

Value 

Units 

Boltzmann  constant 

k 

1.3807  x  10“23 

J  K  - 1 

Elementary  charge 

e 

1.6022  x  10-19 

C 

Electron  mass 

me 

9.1094  x  10-31 

kg 

Proton  mass 

TTlp 

1.6726  x  10-27 

kg 

Gravitational  constant 

G 

6.6726  x  10-11 

m3s-2kg-1 

Planck  constant 

h 

6.6261  x  10“34 

Js 

p* 

II 

?■ 

to 

* 

1.0546  x  10~34 

Js 

Speed  of  light  in  vacuum 

c 

2.9979  x  10“ 

„  1 
m  s 

Permittivity  of 

( 0 

8.8542  x  10-12 

Fra'1 

free  space 

Permeability  of 

Ho 

4tt  x  10~7 

H  rn"1 

free  space 

Proton/electron  mass 

mp/me 

1.8362  x  103 

ratio 

Electron  charge/mass 

efme 

1.7588  x  1011 

Ckg-1 

ratio 

Rydberg  constant 

_  4 

me 

8eo2c/i3 

1.0974  x  107 

m  ~ 1 

Bohr  radius 

ao  =  €oh2 /nme2 

5.2918  x  10-11 

m 

Atomic  cross  section 

r  a02 

8.7974  x  10~21 

m2 

Classical  electron  radius 

re  =  e2 /Income2 

2.8179  x  10-15 

m 

Thomson  cross  section 

(87r/3)re2 

6.6525  x  10-29 

m2 

Compton  wavelength  of 

h/mcc 

2.4263  x  10“12 

m 

electron 

h/me  c 

3.8616  x  10-13 

m 

Fine-structure  constant 

a  =  e2 /2cohc 

7.2974  x  10"3 

a"1 

137.04 

First  radiation  constant 

Ci  =  27 xhc2 

3.7418  x  10-2 

Wm2 

Second  radiation 

C2  =  he/ k 

1.4388  x  10-2 

m  K 

constant 

Stefan -Boltzmann 

o 

5.6705  x  10"8 

W  m  ~  2 1<  “  4 

constant 
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Physical  Quantity 

Symbol 

Value 

Units 

Wavelength  associated 

Aq  =  hc/e 

1.2398  x  10~6 

m 

with  1  eV 

Frequency  associated 

o 

II 

a 

p- 

2.4180  x  1014 

Hz 

with  1  eV 

Wave  number  associated 

II 

O 

8.0655  x  105 

in"1 

with  1  eV 

Energy  associated  with 

hu0 

1.6022  x  10~19 

J 

1  eV 

Energy  associated  with 

he 

1.9864  x  10~25 

J 

1  m-1 

Energy  associated  with 

me3 /8«o2  h2 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

k/e 

8.6174  x  10~5 

eV 

1  Kelvin 

Temperature  associated 

efk 

1.1604  x  104 

K 

with  i  eV 

Avogadro  number 

na 

6.0221  x  1023 

mol'1 

Faraday  constant 

F  =  NAe 

9.6485  x  104 

C  mol-1 

Gas  constant 

R  =  N  a  k 

8.3145 

J  K-1mol~1 

Loschmidt’s  number 

n0 

2.6868  x  1025 

m-3 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6605  x  10  “27 

kg 

Standard  temperature 

To 

273.15 

K 

Atmospheric  pressure 

Po  =  nokTo 

1.0133  x  10s 

Pa 

Pressure  of  1  mm  Hg 

1.3332  x  102 

Pa 

(1  torr) 

1  Molar  volume  at  STP 

Vo  =  RTo/po 

2.2414  x  10"2 

m3 

Molar  weight  of  air 

Mair 

2.8971  x  10-2 

kg 

calorie  (cal) 

4.1868 

J 

Gravitational 

9 

9.8067 

ms-2 

acceleration 

PHYSICAL  CONSTANTS  (cgs)7 


Physical  Quantity 

Symbol 

Value 

Units 

Boltzmann  constant 

k 

1.3807  x  10~16 

erg/deg  (K) 

Elementary  charge 

e 

4.8032  x  10"10 

statcoulomb 

Electron  mass 

me 

9.1094  x  10-28 

(statcoul) 

g 

Proton  mass 

mp 

1.6726  x  10'24 

g 

Gravitational  constant 

G 

6.6726  X  10-8 

dyne-cm2/g2 

Planck  constant 

h 

6.6261  x  10-27 

erg-sec 

h  =  h/2ir 

1.0546  x  10“27 

erg-sec 

Speed  of  light  in  vacuum 

c 

2.9979  x  1010 

cm/sec 

Proton/electron  mass 

mp/me 

1.8362  x  103 

ratio 

Electron  charge/mass 

e/mc 

5.2728  x  1017 

statcoul/g 

ratio 

Rydberg  constant 

n  2ir2me4 

*°°  “  cfc3 

1.0974  x  105 

_  - 1 
cm 

Bohr  radius 

a0  =  h2 /me2 

5.2918  x  10-9 

cm 

Atomic  cross  section 

2 

IT  do 

8.7974  x  10-17 

2 

cm 

Classical  electron  radius 

_  2  2 
re  =  e  /me 

2.8179  x  10~13 

cm 

Thomson  cross  section 

(8ir/3)re2 

6.6525  x  10-25 

cm2 

Compton  wavelength  of 

h/mec 

2.4263  x  10"10 

cm 

electron 

h/me  c 

3.8616  x  10-11 

cm 

Fine-structure  constant 

a  =  e2 /he 

7.2974  x  10-3 

First  radiation  constant 

a'1 

ci  =  2irhc 2 

137.04 

3.7418  x  10-5 

erg-cm2/sec 

Second  radiation 

C2  —  he/ k 

1.4388 

cm-deg  (K) 

constant 

Stefan- Boltzmann 

a 

5.6705  x  10~5 

erg/cm2- 

constant 

Wavelength  associated 

-^0 

1.2398  x  10-4 

sec-deg4 

cm 

with  1  eV 
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Physical  Quantity 


Symbol 


Value 


U  nits 


Frequency  associated 

vo 

2.4180  x  1014 

Hz 

with  1  eV 

Wave  number  associated 

ko 

8.0655  x  103 

cm  ~  1 

with  1  eV 

Energy  associated  with 

1.6022  x  10-12 

erg 

1  eV 

Energy  associated  with 

1.9864  x  10-16 

erg 

1  cm-1 

Energy  associated  with 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

8.6174  x  10-5 

eV 

1  deg  Kelvin 

Temperature  associated 

1.1604  x  104 

deg  (K) 

with  1  eV 

Avogadro  number 

na 

6.0221  x  1023 

mol  -  1 

Faraday  constant 

F  =  N  Ae 

2.8925  x  1014 

statcoul/mol 

Gas  constant 

R  =  N  A  k 

8.3145  x  107 

erg/deg-mol 

Loschmidt’s  number 

n0 

2.6868  x  1019 

cm-3 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6605  x  10-24 

g 

Standard  temperature 

To 

273.15 

deg  ( K) 

Atmospheric  pressure 

p  o  —  nokTo 

1.0133  x  106 

dyne/cm2 

Pressure  of  1  mm  Hg 

1.3332  x  103 

dyne/cm2 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTo  /  po 

2.2414  x  104 

cm3 

Molar  weight  of  air 

Mair 

28.971 

g 

calorie  (cal) 

4.1868  x  107 

erg 

Gravitational 

9 

980.67 

cm /sec2 

acceleration 
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FORMULA  CONVERSION8 

Here  <*  =  102  cinin-1,  ft  =  107  erg  J_1,  f0  =  8.854  2  x  10"  12  F  in"  1 , 
/to  =  4 ir  x  10-7  H  m~\  c  =  (fofi0)~l/2  =  2.9979  x  10®  ins-1,  and 
h  =  1.0546  x  10~34Js.  To  derive  a  dimensionally  correct  SI  formula 
from  one  expressed  in  Gaussian  units,  substitute  for  each  quantity  ac¬ 
cording  to  Q  =  fcQ,  where  k  is  the  coefficient  in  the  second  column 
of  the  table  corresponding  to  Q  (overbars  denote  variables  expressed  in 
Gaussian  units).  Thus,  the  formula  do  =  h2 /me2  for  the  Bohr  radius 
becomes  aao  =  (hft)2 /{(mft  /  a2)(e2  otft  /  4  7r  c0 )] ,  or  a  o  =  e0h2  /  nmc2 .  To 
go  from  SI  to  natural  units  in  which  h  =  c  =  1  (distinguished  by  a  cir¬ 
cumflex),  use  Q  =  k~1Q ,  where  k  is  the  coefficient  corresponding  to  Q 
in  the  third  column.  Thus  do  =  4-k e0h2 /{(mh / c)(e2  cohc)]  —  47r /me2. 
(In  transforming  from  SI  units,  do  not  substitute  for  to>  Mo,  or  c.) 


Physical  Quantity 

Gaussian  Units  to  SI 

Natural  Units  to  SI 

Capacitance 

aj  4ir6o 

to-1 

Charge 

(aft/4-ne0)1/2 

(f0ftc)-1/2 

Charge  density 

(P/4-jraS(0)1/2 

(fo  he)-1'2 

Current 

(  aft /4ir  c0) 1/2 

(po/Kc)1'2 

Current  density 

(/3/4‘ir  a3e0)1/2 

(/.o/ftc)1'* 

Electric  field 

{4-Kfttofa3)1^2 

(c0/hc)1/2 

Electric  potential 

(47r/3e0/a)1/2 

(t„/Rc)V2 

Electric  conductivity 

(4ire0)_1 

Co-1 

Energy 

P 

(dc)-1 

Energy  density 

P/a3 

(he)-' 

Force 

P/a 

(he)-1 

Frequency 

1 

c-1 

Inductance 

4  7r  e0  /  a 

Mo  1 

Length 

a 

1 

Magnetic  induction 

{4*p/a3  no)1/2 

(Mod.c)~1/2 

Magnetic  intensity 

(47r/io  P/a3)l/2 

U‘o/hc)'/2 

Mass 

P/a 2 

c/h 

Momentum 

P/a 

h~ 1 

Power 

P 

(he2)-1 

Pressure 

P/a3 

(he)~l 

Resistance 

4ireo/a 

(fo/ Mo)17" 

Time 

1 

c 

Velocity 

a 

c-1 

18 


MAXWELL’S  EQUATIONS 


Name  or  Description 

SI 

Gaussian 

Faraday’s  law 

dB 

V  x  E  = 

dt 

„  „  1  dB 

c  dt 

Ampere’s  law 

dB 

V  x  H  =  +  J 

dt 

„  „  1  3D  4tt 

V  x  H  =  — - - h  — J 

c  dt  c 

Poisson  equation 

V  •  D  =  p 

V  •  D  =  47T/9 

(Absence  of  magnetic 
monopoles] 

V  -  B  =  0 

V  •  B  =  0 

Lorentz  force  on 

q  (E  +  v  x  B) 

/  1  \ 
q  (  E  -f-  —v  X  B  ] 

charge  q 

V  c  / 

Constitutive 

D  =  eE 

D  =  eE 

relations 

B  =  pH 

B  =  pH 

In  a  plasma,  p  «  p0  =  4ir  x  10  7  H  m  1  (Gaussian  units:  ft  «  1). 
The  permittivity  satisfies  e  «  €0  =  8.8542  x  10-12  F  m-1  (Gaussian: 
c  «  I)  provided  that  all  charge  is  regarded  as  free.  Using  the  drift 
approximation  vj_  =  E  x  B f  B2  to  calculate  polarization  charge  density 
gives  rise  to  a  dielectric  constant  K  =  e/co  =  1  +  367r  x  109p/J32  (SI) 
=  1  -f  4tt pc2  /  B2  (Gaussian),  where  p  is  the  mass  density. 

The  electromagnetic  energy  in  volume  V  is  given  by 


=  i  /  *V(H  . 
J  V 

=  h  [ dv(H 

J  v 


B  +  E  D) 

-  B  +  E  •  D) 


(SI) 


(Gaussian). 


Poynting’s  theorem  is 


dW 


dt 


+  /  N  •  dS  =  -  /  dVJ  •  E, 


v 


where  S  is  the  closed  surface  bounding  V  and  the  Poynting  vector  (energy 
flux  across  S)  is  given  by  N  =  E  x  H  (SI)  or  N  =  cE  x  H/47T  (Gaussian). 
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ELECTRICITY  AND  MAGNETISM 


In  the  following,  e  =  dielectric  permittivity,  /x  =  permeability  of 
conductor,  \i'  =  permeability  of  surrounding  medium,  a  =  conductivity, 
/  =  u>/2*  =  radiation  frequency,  =  /x//xo  and  /ce  =  Where  sub¬ 

scripts  are  used,  ‘1*  denotes  a  conducting  medium  and  ‘2*  '  propagating 
(lossless  dielectric)  medium.  All  units  are  SI  unless  otherw.  •*  specified. 


Permittivity  of  free  space 
Permeability  of  free  space 


Resistance  of  free  space 

Capacity  of  parallel  plates  of  area 
A,  separated  by  distance  d 

Capacity  of  concentric  cylinders 
of  length  /,  radii  a,  b 

Capacity  of  concentric  spheres  of 
radii  a,  b 

Self-inductance  of  wire  of  length 
l ,  carrying  uniform  current 

Mutual  inductance  of  parallel  wires 
of  length  l,  radius  a,  separated 
by  distance  d 

Inductance  of  circular  loop  of  radius 
6,  made  of  wire  of  radius  a, 
carrying  uniform  current 

Relaxation  time  in  a  lossy  medium 

Skin  depth  in  a  lossy  medium 

Wave  impedance  in  a  lossy  medium 

Transmission  coefficient  at 
conducting  surface9 
(good  only  for  T  <  1) 

Field  at  distance  r  from  straight  wire 
carrying  current  I  (amperes) 

Field  at  distance  z  along  axis  from 
circular  loop  of  radius  a 
carrying  current  I 


€0  =  8.8542  x  10-12  Fm-1 

Ho  =  4-k  x  10~~7  H  m ~ 1 
=  1.2566  x  10-6  Hm_1 

R0  =  (no/€0)1/2  =  376.73  ft 

C  =  eA/d 

C  =  2T€/ln(6/a) 

C  =  4xea6/(6  —  a) 

L  =  fil 

L  =  (fi'l/4-K)  [1  +  4  ln(d/a)j 

L  —  b  {/x'  [ln(86/a)  —  2]  +-  ft /4^ 
r  =  (/a 

6  =  (2/ujfia)1^2  =  (x/^<r)~1/2 

Z  =  +  i<r/u>)]1/2 

T  =  4.22  x  10-4(fKTnlKe2/a)1/2 

Be  —  pi f2*T  tesla 

=  0.2//rgauss  (r  in  cm) 

Bz  =  fia2 1 /[2(a2  +  z2f'2) 
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ELECTROMAGNETIC  FREQUENCY/ 
WAVELENGTH  BANDS10 


Designation 

Frequency  Range 

Wavelength  Range 

Lower 

Upper 

Lower 

Upper 

ULF* 

10  Hz 

3  Mm 

ELF* 

10  Hz 

3  kHz 

100  km 

VLF 

3  kHz 

30  kHz 

10  km 

LF 

30  kHz 

300  kHz 

1  km 

MF 

300  kHz 

3  MHz 

100  m 

HF 

3  MHz 

30  MHz 

10  m 

VHF 

30  MHz 

300  MHz 

1  m 

10  m 

UHF 

300  MHz 

3  GHz 

10  cm 

1  m 

SHFf 

3  GHz 

30  GHz 

1  cm 

10  cm 

S 

2.6 

3.95 

7.6 

11.5 

G 

3.95 

5.85 

5.1 

7.6 

J 

5.3 

8.2 

3.7 

5.7 

H 

7.05 

10.0 

3.0 

4.25 

X 

8.2 

12.4 

2.4 

3.7 

M 

10.0 

15.0 

2.0 

3.0 

P 

12.4 

18.0 

1.67 

2.4 

K 

18.0 

26.5 

1.1 

1.67 

R 

26.5 

40.0 

0.75 

1.1 

EHF 

30  GHz 

300  GHz 

1  mm 

1  cm 

Submillimeter 

300  GHz 

3  THz 

100  /im 

1  mm 

Infrared 

3  THz 

430  THz 

700  nm 

100  /rm 

Visible 

430  THz 

750  THz 

400  nm 

700  nm 

Ultraviolet 

750  THz 

30  PHz 

10  nm 

400  nm 

X  Ray 

30  PHz 

3  EHz 

100  pm 

10  nm 

Gamma  Ray 

3  EHz 

100  pm 

Note:  In  spectroscopy  the  angstrom  (A)  is  sometimes  used  (1  A  =  10  8  cm 
=  0.1  nm). 

*The  boundary  between  ULF  and  ELF  is  variously  defined. 

|The  SHF  (microwave)  band  is  further  subdivided  approximately  as  shown.1 1 
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AC  CIRCUITS 


For  a  resistance  /?,  inductance  L,  and  capacitance  C  in  series  with 
a  voltage  source  V  =  Vo  exp(tu/t)  (here  t  =  y/ —  1),  the  current  is  given 
by  I  =  dq/dt,  where  q  satisfies 


L 


=  V. 


Solutions  are  q(t)  =  qa  +-  qt ,  I(t)  =  I,  +  Ity  where  the  steady  state  is 
Ia  =  iu>qa  =  V/Z  in  terms  of  the  impedance  Z  =  R  - f  i(u>L  —  1  /a»C)  and 
/t  =  dqt/dt.  For  initial  conditions  ^(0)  =  90  =  ?o  +  /(0)  =  Jo,  the 

transients  can  be  of  three  types,  depending  on  A  =  R 2  —  4 L/C: 

(a)  Overdamped,  A  >  0 


Qt 

It 


Jo +  7+tfo  ,  ^  /o +7-4o  , 

exp(  — 7_  t) - exp(-7+t), 


7+  “  7- 
7+(/o  +  7-4o) 
7+  “  7- 


exp(-7+t)  - 


7+  ~  7- 

7- (Jo  +  7+4o) 


7+  ~  7- 


exp(-7_t), 


where  7+  =  (R  rfc  AxJ2)/2L; 

(b)  Critically  damped,  A  =  0 


qt  -  [40  +  (Jo  +  7K4o)t]  exp(-7*t), 

It  =  [Jo  -  (Jo  +  7R4o)7«tl  exp(-7/?t), 


where  7n  =  R/2L; 

(c)  Underdamped,  A  <  0 


4t  = 

It  = 


exp(-7i?t), 


7R4o  +  Jo  .  .  ,  _ 

- sin  u\  t  +-  qo  cos  u>i  t 

wi 

'  .  ("I2  +  lR2)qo  +  IrIo  .  . 

io  coswit  -  - sin(u>it) 

Wi 


exp(  —  7/?t), 


Here  wi  =  u/o(l  —  R2C/4L)XI2,  where  u>o  =  ( LC)~X^ 2  is  the  resonant 
frequency.  At  u)  =  u>q,  Z  —  R.  The  quality  of  the  circuit  is  Q  =  u>qL/R. 
Instability  results  when  L,  H,  C  are  not  all  of  the  same  sign. 
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IMENSIONLESS  NUMBERS  OF  FLUID  MECHANICS12 


Name(s) 

Symbol 

Definition 

Significance 

Alfven, 

Karman 

Al,  Ka 

vA/v 

^(Magnetic  force/ 
inertial  force)1/2 

Bond 

Bd 

(. P '  ~P)L29/Z 

Gravitational  force/ 
surface  tension 

Boussinesq 

B 

V/(2  gR)'/7 

(Inertial  force/ 

gravitational  force)1/2 

Brinkman 

Br 

pV2/kAT 

Viscous  heat/conducted  heat 

Capillary 

Cp 

nv/z 

Viscous  force/surface  tension 

Carnot 

Ca 

(T2  -  Tx)/T* 

Theoretical  Carnot  cycle 
efficiency 

Cauchy, 

Hooke 

Cy,  Hk 

pv2/r  =  m2 

Inertial  force/ 

compressibility  force 

Chandra¬ 

sekhar 

Ch 

B2  L2  / pvt] 

Magnetic  force/dissipative 
forces 

Clausius 

Cl 

LV3p/kAT 

Kinetic  energy  flow  rate/heat 
conduction  rate 

Cowling 

C 

(VA/V)2  =  Al2 

Magnetic  force/inertial  force 

Crispation 

Cr 

/x/c/E  L 

Effect  of  diffusion/effect  of 
surface  tension 

Dean 

D 

D3/2V7i/(2r)1/2 

Transverse  flow  due  to 

curvature/longitudinal  flow 

(Drag 

coefficient) 

Cd 

( P '  ~  P)l9 / 
pV2 

Drag  force/inertial  force 

Eckert 

E 

V2  /  cpAT 

Kinetic  energy/change  in 
thermal  energy 

Ekman 

Ek 

(u/2  DL2)1/2  = 
(Ro/Re)l/2 

(Viscous  force/Coriolis  force)1/2 

Euler 

Eu 

A  p/pV2 

Pressure  drop  due  to  friction/ 
dynamic  pressure 

Froude 

Fr 

V/(sL)1'2 

V/NL 

t(Inertial  force/gravitational  or 

1  /  O 

buoyancy  force)  ' 

Gay-Lussac 

Ga 

1/PAT 

Inverse  of  relative  change  in 
volume  during  heating 

Grashof 

Gr 

9L3PAT/u 2 

Buoyancy  force/viscous  force 

[Hall 

coefficient] 

Ch 

A/ rL 

Gy  rofrequency/ 

collision  frequency 

*(t)  Also  defined  as  the  inverse  (square)  of  the  quantity  shown. 
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Name(s) 


Symbol 


Definition 


Significance 


Hartmann 

H 

BL /(/,,)*/*  = 
(Rm  ReC)1/2 

(Magnetic  force/ 

dissipative  force)1/2 

Knudsen 

Kn 

X/L 

Hydrodynamic  time/ 
collision  time 

Lorentz 

Lo 

Vic 

Magnitude  of  relativistic  effects 

Lundquist 

Lu 

HoLVa/v  = 

A1  Rm 

J  x  B  force/resistive  magnetic 
diffusion  force 

Mach 

M 

V/Cs 

Magnitude  of  compressibility 
effects 

Magnetic 

Mach 

Mm 

V/VA  =  Al-1 

(Inertial  force/magnetic  force)1/2 

Magnetic 

Reynolds 

Rm 

Ho  LV/rj 

Flow  velocity/magnetic  diffusion 
velocity 

Newton 

Nt 

F/pL2  V2 

Imposed  force/inertial  force 

Nusselt 

N 

aL/k 

Total  heat  transfer/thermal 
conduction 

Peclet 

Pe 

LV/k 

Heat  convection/heat  conduction 

Poisseuille 

Po 

D2Ap/pLV 

Pressure  force/viscous  force 

Prandtl 

Pr 

v/k 

Momentum  diffusion/ 
heat  diffusion 

Rayleigh 

Ra 

gH3pAT/vK 

Buoyancy  force/diffusion  force 

Reynolds 

Re 

LV/v 

Inertial  force/ viscous  force 

Richardson 

Ri 

( NH/AV )2 

Buoyancy  effects/ 

vertical  shear  effects 

Rossby 

Ro 

V/2QL  sin  A 

Inertial  force/Coriolis  force 

Stanton 

St 

al  pCpV 

Thermal  conduction  loss/ 
heat  capacity 

Stefan 

Sf 

cLT3 /k 

Radiated  heat/conducted  heat 

Stokes 

s 

u!L2f 

Viscous  damping  rate/ 
vibration  frequency 

Strouhal 

Sr 

SL/V 

Vibration  speed/flow  velocity 

Taylor 

Ta 

(2  QL2/u)2 
R1/2(AR)3/2 
(^/* ') 

Centrifugal  force/viscous  force 
(Centrifugal  force/ 
viscous  force)1/2 

Thring, 

Boltzmann 

Th,  Bo 

pcp  V /  ecrT3 

Convective  heat  transport/ 
radiative  heat  transport 

Weber 

W 

pLV2 /E 

Inertial  force/surface  tension 

lomenclature: 

Magnetic  induction 

c 

Speeds  of  sound,  light 

Specific  heat  at  constant  pressure  (units  m2s-2K-1) 

-  2  R 

Pipe  diameter 

Imposed  force 

Vibration  frequency 

Gravitational  acceleration 

,  L 

Vertical,  horizontal  length  scales 

=  PCP  K 

Thermal  conductivity  (units  kgm-1  s~2) 

r  =  (g/H)1/2 

Brunt-Vaisala  frequency 

Radius  of  pipe  or  channel 

Radius  of  curvature  of  pipe  or  channel 

L 

Larmor  radius 

1 

Temperature 

Characteristic  flow  velocity 

A  =  B/W)I/2 

Alfven  speed 

dT 

Newton’s-law  heat  coefficient,  k— —  =  qAT 

dx 

Volumetric  expansion  coefficient,  dV/V  =  (3dT 

Bulk  modulus  (units  kg  m-1  s-2) 

iR,  AV,  Ap,  AT 

Imposed  difference  in  two  radii,  velocities, 
pressures,  or  temperatures 

Surface  emissivity 

Electrical  resistivity 

Thermal  diffusivity  (units  m2s_1) 

Latitude  of  point  o:i  earth’s  surface 

Collisional  mean  free  path 

=  P^ 

Viscosity 

0 

Permeability  of  free  space 

Kinematic  viscosity  (units  m2s_1) 

Mass  density  of  fluid  medium 

/ 

Mass  density  of  bubble,  droplet,  or  moving  object 

Surface  tension  (units  kgs-2) 

Stefan-  Boltzmann  constant 

i 

Solid-body  rotational  angular  velocity 
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SHOCKS 


At  a  shock  front  propagating  in  a  magnetized  fluid  at  an  an^le  6 
with  respect  to  the  magnetic  induction  B,  the  jump  conditions  are  1  ' 

(1)  pU  =  pU  =  q ; 

(2)  pU 2  +  p  +  B  2  /2p  —  pU 2  -f-  p  4-  B j_2  /2/z; 

(3)  pUV  -  £||  B±  /  p.  =  pUV  -  B„Bx/p; 

(4)  B„  =  % 

(5)  BBj.  -  VB|j  =  C/Bx  -  V B„ ; 

(6)  ±(£/2  +  V2)  +  W  +  ( U Bx2  -  yBHBx)/ppB 

=  +  V2)  +  u>  +  (t/Bj2  -  VB^B±)/p.pU. 

Here  U  and  V  are  components  of  the  fluid  velocity  normal  and  tangential 
to  the  front  in  the  shock  frame;  p  =  1/v  is  the  mass  density;  p  is  the 
pressure;  Bj_  =  B  sin  0,  Bj|  =  Bcosfl;  p  is  the  magnetic  permeability 
(p  =  4ir  in  cgs  units);  and  the  specific  enthalpy  is  w  =  e  -f  pt»,  where 
the  specific  internal  energy  e  satisfies  de  =  T ds  —  pdv  in  terms  of  the 
temperature  T  and  the  specific  entropy  s.  Quantities  in  the  region  behind 
(downstream  from)  the  front  are  distinguished  by  a  bar.  If  B  =  0,  then15 

(7)  u  -U  =  [(P-T\v  -v))^2; 

(8)  (p  -  p)(v  -  v)~1  =  q2; 

(9)  w  -  w  =  ±{p  -  p)(v  v), 

(10)  e  -  e  =  |(p  +  p)( v  -  v). 

In  what  follows  we  assume  that  the  fluid  is  a  perfect  gas  with  adiabatic 
index  7  =  1  -f  2/n,  where  n  is  the  number  of  degrees  of  freedom.  Then 
p  =  pRT / m ,  where  R  is  the  universal  gas  constant  and  m  is  the  m^iar 
weight;  the  sound  speed  is  given  by  C,2  =  ( dp/dp)3  =  ~fpv\  and  w  — 
7e  =  7pt»/(7  T  1).  For  a  general  oblique  shock  in  a  perfect  gas  the 
quantity  X  =  t~1{U/Va)2  satisfies14 

(11)  {X  —  ft / a)(X  —  cos2  ff)2  =  X  sin2  8  {[1  +  (r  -  l)/2a]  X  -  cos2  fl}, 
where  r  =  p / p,  ac  =  ^  [7  +  1  —  (7  —  l)r],  and  (3  =  C\,  2  /  V.\  2  =  4  777 p  /  B2 . 
The  density  ratio  is  bounded  by 

(12)  1  <  r  <  (7  +  l)/(7  -  1). 

If  the  shock  is  normal  to  B  (i.e.,  if  6  ~  7r/2),  then 

(13)  U 2  =  (r/ct)  {CB2  +  V a  2  [  1  +  (1  -  7  /  2 )  ( r  -  1)]}; 
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(14)  U/U  =  B/B  =  r; 

(15)  V  =  V; 


(16)  p  -  p  +  (1  -  r~x)pU2  +  (1  -  r2)B2  /2/a. 

If  9  =  0,  there  are  two  possibilities:  switcli-on  shocks,  which  require 
/?  <  1  and  for  which 

(17)  U2  =  rVA2; 

(18)  0  =  Va2/U ; 

(19)  Bx2  =  2B„2(r  -  l)(a  - /?); 


(20)  V  = 

(21)  p  =  p  +  pU2(  1  -  a  +  /?)(1  -  r-1), 

and  acoustic  (hydrodynamic)  shocks,  for  which 

(22)  V 2  =  (r/a)C.2; 

(23)  &  =  (7/rj 

(24)  V=  B±=  0; 

(25)  p  =  p  +  pC/2(l  -T-1). 

For  acoustic  shocks  the  specific  volume  and  pressure  are  related  by 

(26)  v/v  =  [(7  +  l)p  -f  (7  -  l)p]  /  ((7  -  l)p  +  (7  +  1)p]- 
In  terms  of  the  upstream  Mach  number  M  =  U/Ce, 

(27)  p/p  =  v/v  =  U/U  =  (7  +  l)M2/[(7  -  1)M2  +  2]; 

(28)  p/p  =  (27M2  7+l)/(7+l); 

(29)  T/T  =  [(7  -  1  )M2  +  2](27Af 2  -  7  +  l)/(7  +  1)2M2; 

(30)  M2  =  [(7  -  1)M2  +  2]/[27M2  -7  +  1]. 


The  entropy  change  across  the  shock  is 

(31)  As  =  3  -  s  =  cv  In [(p/p)(p/py]. 


where  cv  —  R/( 7  —  l)m  is  the  specific  heat  at  constant  volume;  here  R 
is  the  gas  constant.  In  the  weak-shock  limit  (M  — »  1), 


/oox  a  27(7  -  1)  ,  ,r2  ,  \ 3 

(32)  As  —  cv  — - ~~(M  -1) 

3(7  +  1) 


I67  R 


3(7  +  l)m 


(M  -  1)' 


The  radius  at  time  t  of  a  strong  spherical  blast  wave  resulting  from  the 
explosive  release  of  energy  E  in  a  medium  with  uniform  density  p  is 


(33)  Rs  -  Co(Et2/p)x>\ 


where  Co  is  a  constant  depending  on  7.  For  7  =  7/5,  Co  =  1.033. 
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FUNDAMENTAL  PLASMA  PARAMETERS 


All  quantities  are  in  Gaussian  cgs  units  except  temperature  (T,  Te, 
Tj)  expressed  in  eV  and  ion  mass  (m,)  expressed  in  units  of  the  proton 
mass,  ii  —  m.ijmp ;  Z  is  charge  state;  k  is  Boltzmann’s  constant;  K  is 
wavelength;  7  is  the  adiabatic  index;  In  A  is  the  Coulomb  logarithm. 


Frequencies 

electron  gyrofrequency 

ion  gyrofrequency 


electron  plasma  frequency 


ion  plasma  frequency 


electron  trapping  rate 

ion  trapping  rate 

electron  collision  rate 
ion  collision  rate 


/„  =  u>ce/2it  =  2.80  x  106B  Hz 

utce  =  eB /me  c  =  1.76  x  107  B  rad /sec 

fd  =  u>ci/2x  =  1.52  x  103Zn~1B  Hz 

u >ci  =  Z eB / TTLiC  —  9.58  X  10 3  Z fi~1  B  rad /sec 

fpe  =  u>pe/2-K  =  8.98  X  103ne1/2  Hz 

u>pc  =  (4?r  nee2/me)1/2 

=  5.64  x  104ne  1^2  rad/sec 

f pi  =  Wpi /2tt 

=  2.10  X  10 2ZfjL~ll2ni1/2  Hz 
ujpi  —  (AtruiZ2  e2 /rm)1^2 

=  1.32  X  103Z/i"1/2nt1/2rad/sec 
t'Te  =  (eKEfme)1/2 

=  7.26  X  10 bK1/2E1/2  i  ec  1 
vTi  =  ( ZeKE/mi )1/2 

=  1.69  X  107 Z1/2  K1/2  E1/2 n~1/2  sec-1 
ve  =  2.91  x  10_6ne  In  A Te~3/2  sec-1 
Vi  =  4.80  x  10_8Z4/i_1/2n,  In  A Ti~3/2  sec-1 


Lengths 

electron  deBroglie  length 

classical  distance  of 
minimum  approach 

electron  gyroradius 
ion  gyroradius 

plasma  skin  depth 
Debye  length 


X  -  h/(mckTe)1/2  =  2.76  x  10-8Te~1/2  cm 
e2 / kT  =  1.44  x  10-7T_1  cm 

re  =  vj'c  f  u>ce  =  2.38Te1^2  B  1  cm 

Ti  =  VTi/Ud 

=  1.02  x  10 2/i1/2Z~1T,1/2B"1  cm 
c / u>pC  —  5.31  X  10snc-1^2  cm 
XD  =  (kT/4*ne2)1/2 

=  7.43  x  102r1/2n-1  /2  cm 
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Velocities 

electron  thermal  velocity  t>Te  =  (kTe  /me)1^2 

=  4.19  x  10 7T*e1/2  cm/scc 

ion  thermal  velocity  vri  =  (kTi  /rrii)1^2 

=  9.79  x  105/i“  1/2 Ti 1/2  cm/scc 
ion  sound  velocity  C*  =  (7 Z  kTe /m,i)1^2 

=  9.79  X  105(7ZTc//i)1/2  cm  /see 
Alfven  velocity  v ^  —  B/(47rn,  mj 

=  2.18  x  101V_1/2ni_1/2Bcm/sec 

Dimensionless 

(electron/proton  mass  ratio)1/2  (me  jmv)x^2  =  2.33  X  10-2  =  1/42.9 

number  of  particles  in  (47r/3)nA£> 3  =  1.72  x  109T3/2n~ 1/2 

Debye  sphere 

Alfven  velocity/speed  of  light  v a  / c  =  7.2 8/i-1/2n;  -1//2B 

electron  plasma/gyrofrequency  vpt./uice  =  3.21  x  10~3nc  1/2  B"1 
ratio 

ion  plasma/gyrofrequency  ratio  upi/vd  =0.137 /iX/2n, 1//2 B-1 
thermal/magnetic  energy  ratio  (5  —  SrnkT  /  B2  —  4.03  x  10-11nTB~2 
magnetic/ion  rest  energy  ratio  B2/87rn;m,c2  =  26.5^t -1n, -1B2 

Miscellaneous 

Bohm  diffusion  coefficient  Dg  =  ( ckT/16eB ) 

=  6.25  X  106TB-1  cm2/sec 

transverse  Spitzer  resistivity  77^  =  1.15  x  10-14Z  In  AT-3/2  sec 

=  1.03  x  10~2Zln  AT~3/2  Qcm 

’he  anomalous  collision  rate  due  to  low-frequency  ion-sound  turbulence 
> 

v*  «  cupeW/kT  =  5.64  X  104ne  1/2 W" / fcT  sec  1 , 

/here  W  is  the  total  energy  of  waves  with  u> / K  <  vxi- 
lagnetic  pressure  is  given  by 

Pmag  =  B2/87r  =  3.98  x  106B2  dynes/cin2  =  3.93(B/Bq)2  atm, 

/here  Bo  =  10  kG  =  1  T. 

)etonation  energy  of  1  kiloton  of  high  explosive  is 

Wkx  —  1012  cal  =  4.2  x  1019  erg. 
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PLASMA  DISPERSION  FUNCTION 
Definition16  (first  form  valid  only  for  Im(  >  0): 

'  +  0°  dt  exp  (~f2) 


2(0 


t  -  c 


=  2*  exp 


(-<2)  f  dt  exp  (-<2)  • 

v  —  OO 


Physically,  £  =  x  +  iy  is  the  ratio  of  wave  phase  velocity  to  thermal 
velocity. 


Differential  equation: 
dZ 


d< 


=  -2(1  +  CZ)  ,  Z( 0)  =  iTr 


-  ,'1/2. 


d2Z  dZ 

— -  +  2< —  +  2Z  =  0. 

dC  d( 


Real  argument  ( y  =  0): 

Z(x)  =  exp  (-z2)  ^i'Jr1/^2  —  2  J  dt  exp  (f 

Imaginary  argument  (x  =  0): 

Z(iy )  =  »*1/2  exp  (y2)  [1  -  erf(y)]  . 

Power  series  (small  argument): 

2(0  =  *^1/2  exp  (-C2)  -  2<  (l  -  2<2/3  +  4(4/15  -  8C6/105  +  ••-). 
Asymptotic  series,  |£|  1  (Ref.  17): 

2(0  ~  ^1/2<texp  (-C2)  -  C"1  (l  +  1/2C2  +  3/4C4  +  15/8C6  +  •••), 
where 


~  i 


y  >  1*1 
ivi  <  m-1 


1-1 


y  <  -|*| 

Symmetry  properties  (the  asterisk  denotes  complex  conjugation): 

z«*)  =  -(Z(-<)]*; 

Z«*)  =  [Z(OJ  *  +  2iir1/2  exp[-«*)2l  (j,  >  0). 

Two-pole  approximations18  (good  for  (  in  upper  half  plane  except  when 

y  <  7T1/2  x2  exp(  —  x2),  x  1): 

,  0.50  +  0.81*  0.50  —  0.81i 

Z(Q  «  - 7 - — — ,  a  =  0.51  -  0.81*; 


a  -  C 


a*  +  C 


,,  0.50  +  0.96*  0.50  -  0.96* 

Z  (C)  « -  +  — — — - — ,  b  =  0.48  —  0.91*. 


(6-C)2 


(&*  +  C)2 
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COLLISIONS  AND  TRANSPORT 


Temperatures  are  in  eV;  the  corresponding  value  of  Boltzmann’s  constant 
fc  =  1.G0  X  10~12erg/eV;  masses  /r,  0.'  are  in  units  of  the  proton  mass; 
v  =  Za  e  is  the  charge  of  species  c*.  All  other  units  are  cgs  except  where 
oted. 

telaxation  Rates 

Rates  are  associated  with  four  relaxation  processes  arising  from  the  in¬ 
fraction  of  test  particles  (labeled  a)  streaming  with  velocity  va  through  a 
ackground  of  field  particles  (labeled  /3): 

slowing  down 

a  _  _.a/0__ 

dt  ~  U*  Vq 

transverse  diffusion 

^  /  -  \2  a/0  2 

—  (v„  -  v„)x  =  tt±  va 

parallel  diffusion 

d  /  -  \2  a/0  2 

*(v»  -v«)||  ="/  «• 

energy  loss 

^  2  a/0  2 

diV°  =  ’ 

vhere  the  averages  are  performed  over  an  ensemble  of  test  particles  and  a 
Maxwellian  field  particle  distribution.  The  exact  formulas  may  be  written19 


"s10  =  (1  +  rnQ /mt3)iP(xa/0)u°/0\ 

»l/0  =  2  [(1  -  l/2xa/0)iP(xa/0)  +  rp'(xQ/0)]  Vq/0 \ 

I'll  =  [*{x  ,P)/x  /PJ  vQ  ; 

v*/0  =  2  [(m0  /m0)Tp(xn/0)  -  ip'{xa/0)^  v°/0 , 


vhere 


u 


a  /  0 
0 


2  2 

A-kco,  eft  A^jnj/m, 


2  3 


X,y/ 0  ~  mf1va2 /2kT;1; 


V’(x)  = 


{x) 


dip 

dx 


tml  A,k  i  —  Id  A„,<  is  the  Couloml)  logarithm  (see  below).  Limiting  forms  of 
,  j/j_  and  i'||  art?  given  in  the  following  table.  All  the  expressions  shown 


31 


have  units  cm3  sec-1.  Test  particle  energy  e  and  field  particle  temperature  T 
are  both  in  eV;  /x  =  mi/mp  where  mp  is  the  proton  mass;  Z  is  ion  charge 
state;  in  electron-electron  and  ion-ion  encounters,  field  particle  quantities  are 
distinguished  by  a  prime.  The  two  expressions  given  below  for  each  rate  hold 
for  very  slow  (xa^0  <£  1)  and  very  fast  (xa^0  1)  test  particles,  respectively. 


Electron-electron 


v\  /ne/Aee/  «  5.8  x  10  T 


Slow 


—  6  rp  —  3  /  2 


Fast 


^l/e  /ne/Aee/  «  5.8  x  10-6T-1/2e-1  — ♦  7.7  x  10-o€-J/z 

v*/e>  In.t  A.,/  «  2.9  x  10-6T-1/2c-1  — ♦  3.9  x  10-6Te-5/2 

II  e  ee 

Electron-ion 

u]lifniZ2\ei  «  0.23 /i3/2T-3/2  — ►  3.9  x  10"6e'3/2 

/niZ2\ci  «  2.5  x  10-4/x1/2T-1/2e-1 — *  7.7  x  10-6e-3/2 
v*.fi  /mZ2\ei  »  1.2  x  10-4/i1/2T-1/26-1 — ♦  2.1  x  10-9/x-1Te-5/2 


7.7  x  10-6f-3/2 
7.7  x  10-6e“3/2 
3.9  x  10-6Te-5/2 


3.9  x  10-6c-3/2 


Ion-electron 

t/*/e fneZ2Xie  «  1.6  x  10-9/x_1T-3/2 


i S/€/neZ2\u  «  3.2  x  10-5V_1T-1/*e-1 
v\/*lntZ2  Aie  «  1.6  x  10-9/x_1T-1^2e-1 


-9  u-1T-1/2e-1 


1.7  x  10-4/i1/2e-3/2 

1.8X  10-V1/2«"3/2 

1.7  x  10-4/x1/2Te-5/2 


Ion-ion 


ni/Z2Z/2A.-, 


a.8xl0-.^f1+£l)T 

M  V  ^  / 


-3/2 


9.0  x  10' 


(;♦?) 


A*  ' 

63/2 


n,/Z2Z/2Ait/ 


1  „  w  .  n-r  11/2  -lm-l/2  -1 
«  1.4  x  10  /x  '  /x  T  c 


1.8  x  10-7/x- 1/2 c-3/2 


n,/Z2Z'2  A,,/ 


6.8  x  10 -8/x,1/2/x-1T-1/2e  1 


-.9.0X  10-8/x1/2/-1Te-5/2 

In  the  same  limits,  the  energy  transfer  rate  follows  from  the  identity 

=  2i/s  -  i/j_  -  ^h, 

except  for  the  case  of  fast  electrons  or  fast  ions  scattered  by  ions,  where  the 
leading  terms  cancel.  Then  the  appropriate  forms  are 

vcJ'  - ♦  4.2  x  10-9n,Z2Aei 

[<r3/2/x-1  -  8.9  x  104(/x/T)1/2e-1  cxp(  —  1836/xt/T)j  sec-1 
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and 


v\/x'  — *  1.8  x  10_7n,/Z2Z'2AIJ/ 

[«  3/2/i1/2//x/  -  l.l(/i,/T)1/2e-1  exp(-^/f/T)j  see-1. 

In  general,  the  energy  transfer  rate  i/® 1 0  is  positive  for  c  >  ea*  and  nega¬ 
tive  for  t  <  eQ*,  where  x*  =  (tu0 /mQ)ea* /T0  is  the  solution  of  xj>' (x*)  = 
(mo  /m0 )V;(z*).  The  ratio  eQ* /T0  is  given  for  a  number  of  specific  a,  in  the 
following  table: 


a/p 

i/e  e/e,  t/t  e/p  e/D  e/T,  e/He3  e/He4 

fa* 

1.5  0.98  4.8  x  10-3  2.6  x  10-3  1.8  x  10-3  1.4  x  10-3 

When  both  species  are  near  hdaxwellian,  with  T;  ^  Tc,  there  are  just 
two  characteristic  collision  rates.  For  Z  =  1, 

r/e  =  2.9  X  10-6nATe~3/2  sec-1 ; 

Vi  =  4.8  x  10-8nA7\“3/2p“1/2  sec-1. 


Temperature  Isotropization 
Isotropization  is  described  by 


dTx 

dt 


1  Si 

2  dt 


=  -«4(7x  -T„), 


where,  if  A  =  Tx/T||  -  1  >  0, 


v 


a 

T 


2  e^3  A0^)  ._2 

mrr1/2(fcTj|)3/2 


—  3  +  (A  +  3) 


tan  ^A1/2) 
A1/2 


If  A  <  0,  tan  j(Ai/2)/A1/2  is  replaced  by  tanh  1  (  —  A ) 1  /2 /(  —  A ) 1  ^2  .  For 
T’x  ~  T||  =  T , 


i4  =  8.2  X  10-7nAT-3/2  sec-1; 

*4  =  1-9  x  10-8nAZ2p-1/2T-3/2  sec-1. 
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Thermal  Equilibration 

If  the  components  of  a  plasma  have  different  temperatures,  but  no  rela¬ 
tive  drift,  equilibration  is  described  by 


dTa 

dt 


-  Tn), 


0 


where 


r.n  i  0  —  1  o  ,  r,  —  19  (  mo  ^ 0  )  ^  Za  Z 0  7t-  A  0 

v  —  X.o  X  iu  -  sec 

(ma  T(3  -f  m0Ta  )3/2 

For  electrons  and  ions  with  Te  «  T,  ~  T,  this  implies 


-  i 


v*/l  /  ni  ~  v'J' /ne  =  3.2  x  10_9Z2A/>r3/2  — 3  ---1 


cm  sec 


Coulomb  Logarithm 

For  test  particles  of  mass  ma  and  charge  ea  =  Zae  scattering  off  field 
particles  of  mass  m0  and  charge  e0  =  Z0e,  the  Coulomb  logarithm  is  defined 
as  A  =  In  A  =  ln(rmax/rmin).  Here  rmin  is  the  larger  of  eae^/ma()u2  and 
h/2mo0U,  averaged  over  both  particle  velocity  distributions,  where  ma0  = 

mam0 /(m„  -f  m0)  and  u  =  vQ  -  v0-,  rmax  =  (4*  ^  nyey2  /  kTy)~x/2 ,  where 
the  summation  extends  over  all  species  7  for  which  u2  <  vTy 2  =  kTy/my.  If 
this  inequality  cannot  be  satisfied,  or  if  either  uujca~ 1  <  rmax  or  uuic0~ 1  < 
71  ma.x ■>  the  theory  breaks  down.  Typically  A  ^  10-20.  Corrections  to  the  trans¬ 
port  coefficients  are  0(A_1);  hence  the  theory  is  good  only  to  ~  10%  and  fails 
when  A  ~  1. 

The  following  cases  arc  of  particular  interest: 

(a)  Thermal  electron  electron  collisions 


A  r( 
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-3/2^ 

T, 

<  10 

eV 

1 

=  24  - 

In  (n//2TF' 

-1)- 

Tc 

>  10 

eV 

)  Ele< 

:lron 

ion  c 

ollisions 

A,,  = 
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(c)  Mixed  ion-ion  collisions 


A  —i  A  ■/ •  —  23  —  In 


ZZ'(n  +  n') 

+  /^Ti 


( 


n,Z2 

T, 


+ 


n,  /  Z 
T,/ 


i  2 


1/2"! 


(d)  Counterstreaming  ions  (relative  velocity  vd  —  Pdc)  in  the  presence  of 
warm  electrons,  kTx  / mx ,  fcTp  /rnii  <  vq2  <  kTc/me 


A,-,/  =  A,/  ,•  =  35  —  In 


ZZ'(n  4-  ft')  /  ne  \  1/2 

^'/?D2  V7W 


Fokker- Planck  Equation 

Df°  df 


Dt 


dt 


+  v  •  V/  +  F  -  Vv/ 


-(£) 


coll 


where  F  is  an  external  force  field.  The  general  form  of  the  collision  integral  is 

<8/"/ae)„,,  =  -  Y,„  Vv  with 


j“/a  =  it<\<,0c-Zi£- 

m  a 


J  d\'(u2l  —  uu)u  3 

(  — /°,(v)Vv»//Vv')  -  -i-/^(v,)Vv/°(v)) 
L  m0  m0  J 


(Landau  form)  where  u  =  v'  —  v  and  /  is  the  unit  dyad,  or  alternatively, 


Jo/0  =  4x A 


eo2e02 

°0 


|/"(v)Vvff(v)  -  i  Vv  •  [/”(*)VvVvS(v)]  |  , 


where  the  Rosenbluth  potentials  are 


-I 


G(v)  =  v'  )v-d3v' 


//(v)  =  1  + 


rn. 


£  0  I  I  \  ' —  1  1 3  f 

/  (v  )u  a  */  . 
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If  species  a  is  a  weak  beam  (number  and  energy  density  small  compared  with 
background)  streaming  through  a  Maxwellian  plasma,  then 


>/0  _ 


m< 


a  /  0  r  c*  at  /  0  t-t  r 

'B  v/  -  -u..  vv  ■  Vv/ 


mQ  +  mp  ”  '  2  II 

-  (v2l  -  vv)  •  Vv/°. 


B-G-K  Collision  Operator 

For  distribution  functions  with  no  large  gradients  in  velocity  space,  the 
Fokker-Planck  collision  terms  can  be  approximated  according  to 

The  respective  slown.g-down  rates  u^13  given  in  the  Relaxation  Rate  section 
above  can  be  :  for  i/ap,  assuming  slow  ions  and  fast  electrons,  with  e  re¬ 
placed  by  T,  .  For  i/ee  and  t/u,  one  can  equally  well  use  i/_l,  and  the  result 
is  insensitive  to  whether  the  slow-  or  fast-test-particle  limit  is  employed.)  The 
Max welli uns  Fa  and  Fa  are  given  by 


where  n0,  va  and  Ta  are  the  number  density,  mean  drift  velocity,  and  effective 
temperature  obtained  by  taking  moments  of  fQ.  Some  latitude  in  the  definition 
of  Ta  and  va  is  possible;20  one  choice  is  Te  =  Ti,  Ti  —  Tc,  ve  =  v, ,  v,  =  ve. 

Transport  Coefficients 

Transport  equations  for  a  multispecics  plasma: 


dann 

— - - 1-  n„V  ■  =  0; 

dt 


d 


m„  n. 


dt 


=  -vP( 


V 


P  Or  Za  CTi  Q 


1 

E  4-  -  vQ  x  B 

c 


+  Rev ; 
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3  da  kT 

2n«  — +  Po  v  •  v«  =  ~ v  ■  q«»  -  Pa  ;  Vva  +  ga. 

Here  d°  / dt  =  d / dt  +  v Q  ■  V \  pQ  =  nakTa ,  where  fc  is  Boltzmann’s  constant; 
r„  =  and  =  ]T^  Q a 0 1  where  R.a/3  and  Qo,^  are  respectively 

the  momentum  and  energy  gained  by  the  ath  species  through  collisions  with 
the  /3th;  Pa  is  the  stress  tensor;  and  qa  is  the  heat  flow. 

The  transport  coefficients  in  a  simple  two-component  plasma  (electrons 
and  singly  charged  ions)  are  tabulated  below.  Here  ||  and  _L  refer  to  the  di¬ 
rection  of  the  magnetic  field  B  =  b B;  u  =  ve  —  v,  is  the  relative  streaming 
velocity;  ne  =  nt  =  n;  j  =  —  neu  is  the  current;  tvce  =  1.76  x  107Bsec_1  and 
u>ci  =  (me/m,) ujce  are  the  electron  and  ion  gyrofrequencies,  respectively;  and 
the  basic  collisional  times  are  taken  to  be 


r,  = 


3ymT(fcre)3/2  _ . _5^3/2 


4v/2ir  n\e4 


~  3.44  x  10' 


nA 


sec, 


where  A  is  the  Coulomb  logarithm,  and 


r,  = 


3^mr(tn)»/*=a>09xl0,3^  ,/asec_ 


4  yfr-n.  Ae4 


nA 


In  the  limit  of  large  fields  (u>cotrtt  1,  a  =  t,  e)  the  transport  processes  may 
be  summarized  as  follows:2 


momentum  transfer 

R-ei  = 

frictional  force 

Ru  — 

electrical 

O II  = 

conductivities 

■  1 

thermal  force 

R  t  — 

ion  heating 

Qi  = 

electron  heating 

Qe  = 

ion  heat  flux 

q.  = 

ion  thermal 

K.|  = 

conductivities 

II 

electron  heat  flux 

qr  = 

fiictional  heat  flux 

q'u  = 

=  2.0cr  j_  =  2.0 


ne2  re 


m. 


Rr  =  -0.71nV,i(kTe)  - 
3m,  nk 

- —  — (Te  -  T,); 

mi  re 


3n 


2u>ce  Tg 


b  x  V±(kTe); 


q.  =  -^||V||(fcT,)  -  k'±V  ±(kT,)  +  x'Ab  x  Vx(kT,); 
nkTiT,  •  2  nkT{  ■  5  nkT, 


m. 


*J.  = 


m,  uj  frt- 

1  Cl  1 


*A  = 


2m,  u)c , 


3  nkTe 
2  tJ  £  T  f> 


b  x  uj_; 
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thermal  gradient 
heat  flux 


=  -/cf,V|,(ltre)  -  K lvx(fcre)  -  X  Vj _(kTr); 


electron  thermal 
conductivities 


,  nfcTere  nkTc 

/cm  =  3.2 - — :  k'  =4.7 - - 


m. 


meuc^Te 


bnkT 1 


2mru, 


stress  tensor  (either  Pxx  =  —  —  (Wxx  -f  Wvv)  —  ~-{Wxx  —  Wyy)  —  T}^WXV\ 

2  2 


species) 


PWV=  ~^{WIX  +  Wyy)  +  ^(Wrr  -  PVyy)  +  773  ^y; 

z  z 

Pxy=  Pyx  =  +  ^(IV„  -  V^yy); 

z 


P  i:  —  P  rx  —  *l2^iz 

Pyz  =  Pry  =  ~V2^Vz  +  TJ4  Wx  z  ; 

Prr  =  -7?oVVzx 

(here  the  z  axis  is  defined  parallel  to  B); 


ion  viscosity 


:  .  3  nkTi  :  6nkTi 

r] q  =  0.96 nkTiTi;  t?1  =  ~ 2  7?  V2  =  7 —  i 


10u>  2Ti 


V3 


nkTt 

2uci 


i  V4  = 


nkT; 

UJci 


5u>  2T{ 

Cl  1 


electron  viscosity 


.  .  nkTe  nkT, 

Vo  =0.73  nkTere;  vl^O.51——;  V\  =  2.0 


^3 


nkTe  e  nfcTe 

-  o — ;  ^4  = - 

2u/cc  u/ce 

For  both  species  the  rate-of-strain  tensor  is  defined  as 


U>  *■  Te 
ce  c 


dvj  dvk  2 

Wik  -  — -  +  - - 6jkV  •  v. 

dxk  dxj  3 


When  B  =  0  the  following  simplifications  occur: 


Ru  =  ne  j/<T,,;  RT  =  -0.71  nV(fcTe);  q,  =  -kI.V^T,-); 


qfu  =  0. 7lnkTeu\  qeT  =  — /cjj  V(fcTe);  P>fc  =  -  Vo  Wj  k  ■ 

For  u>ceTe  1  a>c,r, ,  the  electrons  obey  the  high-field  expressions  and  the 
ions  obey  the  zero-field  expressions. 

Collisional  transport  theory  is  applicable  when  ( 1 )  macroscopic  time  rates 
of  change  satisfy  d/dt  <C  1/r,  where  r  is  the  longest  collisional  time  scale,  and 
(in  the  absence  of  a  magnetic  field)  (2)  macroscopic  length  scales  L  satisfy  L 
l,  where  l  =  v r  is  the  mean  free  path.  In  a  strong  field,  u/cer  1.  condition 
(2)  is  replaced  by  L(|  >  /  and  L_l  >  x/^V  (P.L  >•  i»  a  uniform  field), 


where  JL||  is  a  macroscopic  scale  parallel  to  the  field  B  and  Lj_  is  the  smaller 
of  D/\V±D\  and  the  transverse  plasma  dimension.  In  addition,  the  standard 
transport  coefficients  are  valid  only  when  (3)  the  Coulomb  logarithm  satisfies 
A  ^  1;  (4)  the  electron  gyroradius  satisfies  re  \p,  or  87T nt.mf.c2  D2 ;  (5) 
relative  drifts  u  =  va  —  vp  between  two  species  are  small  compared  with  the 
thermal  velocities,  i.e.,  u 2  <C  kTa/ma ,  kTp/nip;  and  (6)  anomalous  transport 
processes  owing  to  microinstabilities  are  negligible. 

Weakly  Ionized  Plasmas 

Collision  frequency  for  scattering  of  charged  particles  of  species  a  by 
neutrals  is 

ua  =  n0<j“ /0(kTa/ma  )1/2, 

where  no  is  the  neutral  density  and  cr^°  is  the  cross  section,  typically  ~ 
5  x  10“ 15  cm2  and  weakly  dependent  on  temperature. 

When  the  system  is  small  compared  with  a  Debye  length,  L  < C  A^j,  the 
charged  particle  diffusion  coefficients  are 

Da  —  kTa  j  7TLa  ^a  1 

In  the  opposite  limit,  both  species  diffuse  at  the  ambipolar  rate 

=  -  VcDr  =  (T,  -f  Tg  )£>,  £)„ 

A  /x,  -  /X.  TiDc  +  TcDi  ' 

where  /.ia  =  ea/m0  va  is  the  mobility.  The  conductivity  aa  satisfies  aa  = 

Or  Or  M  Or  • 

In  the  presence  of  a  magnetic  field  B  the  scalars  fx  and  a  become  tensors, 
r=ffa.E  =  afE||  +  «rlEx  +  ^Exb, 

where  b  =  B /D  and 


Or  _  2  f 

&  j}  —  a  J  Tiz (y  v Q  , 

r*  o  2  j  /  2  .  2  \ 

a x  =  a\\  vc.  /(^o  +  ); 

or  ct  1,2.  2  \ 

~  a\\  ^o^ro  /(*'« 

Here  and  a A  are  the  Pedersen  and  Hall  conductivities,  respectively. 
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APPROXIMATE  MAGNITUDES 
IN  SOME  TYPICAL  PLASMAS 


Plasma  Type 

n  cm  3 

T  eV 

ujpc  sec  1 

Aq  cm 

n\o3 

v'ci  sec  1 

Interstellar  gas 

1 

1 

6  x  104 

7  x  102 

— 

4  X  10® 

7  x  10~6 

Gaseous  nebula 

103 

1 

2  x  106 

20 

107 

6  x  10“2 

Solar  Corona 

106 

102 

6  x  107 

7 

4  x  108 

6  x  10“2 

Diffuse  hot  plasma 

1012 

102 

6  x  1010 

7  x  10-3 

4  x  105 

40 

Solar  atmosphere, 

1014 

1 

6  x  1011 

7  x  10_s 

40 

2  x  109 

gas  discharge 

Warm  plasma 

1014 

10 

6  x  1011 

2  x  10-4 

103 

107 

Hot  plasma 

1014 

102 

6  x  1011 

7  x  10~4 

4  x  104 

4  x  106 

Thermonuclear 

1015 

104 

2  x  1012 

2  x  10~3 

107 

5  x  104 

plasma 

Theta  pinch 

1016 

102 

6  x  1012 

7  x  10~5 

4  x  103 

3  x  108 

Dense  hot  plasma 

1018 

102 

6  x  1013 

7  x  10“6 

4  X  102 

2  x  1010 

Laser  Plasma 

1020 

102 

6  x  1014 

7  x  10"7 

40 

2  x  1012 

•  *\*  *  •  ,  22 
The  diagram  (facing)  gives  comparable  information  in  graphical  form. 
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Iog,0  n  (cm'3) 
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-2 


-1 


!og,0  T (eV) 
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IONOSPHERIC  PARAMETERS23 


The  following  tables  give  average  nighttime  values.  Where  two  num¬ 
bers  are  entered,  the  first  refers  to  the  lower  and  the  second  to  the  upper 
portion  of  the  layer. 


Quantity 

E  Region 

F  Region 

Altitude  (km) 

90-160 

160-500 

Number  density  (m~3) 

1.5  x  1010-3.0  x  1010 

5  x  1010-2  x  1011 

Height-integrated  number 
density  (m-2) 

9  x  1014 

4.5  x  1015 

Ion-neutral  collision 
frequency  (s~  1 ) 

2  X  103-102 

0.5-0.05 

Ion  gyro-/collision 
frequency  ratio  Ki 

0.09-2.0 

4.6  x  102-5.0  x  103 

Ion  Pederson  factor 

Ki/(1  +  K,2) 

0.09-0.5 

2.2  X  10_3-2  x  10“4 

Ion  Hall  factor 
*• « 2  /  ( 1  +  k«2) 

8  x  10~4-0.8 

1.0 

Electron-neutral  collision 
.. 

frequency  (s  ) 

1.5  x  104-9.0  x  102 

80-10 

Electron  gyro-/collision 
frequency  ratio 

4.1  x  102-6.9  x  103 

7.8  x  104-6.2  x  105 

Electron  Pedersen  factor 

Ke/(1  +  Ke2) 

2.7  x  10-3-1.5  x  10~4 

10~s-1.5  x  10_c 

Electron  Hall  factor 

XC2/(1  +  Xe2) 

1.0 

1.0 

Mean  molecular  weight 

28-26 

22-16 

Ion  gyrofrequency  (s-1) 

180-190 

230-300 

Neutral  diffusion 

coefficient  (in2s-1) 

30-5  x  103 

105 

The  terrestrial  magnetic  field  in  the  lower  ionosphere  at  equatorial 


lattitudes  is  approximately  Bo  —  0.35  x 
Re  =  6371  km. 


10 


-4 


tesla.  The  earth’s  radius  is 
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SOLAR  PHYSICS  PARAMETERS24 


Parameter 

Symbol 

Value 

U  nits 

Total  mass 

Mq 

1.99  X  1033 

g 

Radius 

Rq 

6.96  x  1010 

cm 

Surface  gravity 

9Q 

2.74  x  104 

cm  s  _  2 

Escape  speed 

V  oo 

6.18  x  107 

cm  s- 1 

Upward  mass  flux  in  spicules 

— 

1.6  x  10~9 

gem-2  s-1 

Vertically  integrated  atmospheric  density 

— 

4.28 

gem-2 

Sunspot  magnetic  field  strength 

flmix 

2500-3500 

G 

Surface  effective  temperature 

T0 

5770 

K 

Radiant  power 

Tq 

3.83  x  1033 

ergs-1 

Radiant  flux  density 

T 

6.28  x  1010 

erg  cm~2s~ 1 

Optical  depth  at  500  nm,  measured 

0.99 

— 

from  photosphere 

Astronomical  unit  (radius  of  earth’s  orbit) 

AU 

1.50  x  1013 

cm 

Solar  constant  (intensity  at  1  AU) 

/ 

1.36  x  106 

erg  cm-2  s-1 

Chromosphere  and  Corona25 


Parameter  (Units) 

Quiet 

Sun 

Coronal 

Hole 

Active 

Region 

Chromospheric  radiation  losses 
(ergem-2  s-1) 

Low  chromosphere 

2  x  106 

2  X  106 

>  io7 

Middle  chromosphere 

2  x  106 

2  X  106 

107 

Upper  chromosphere 

3  x  10s 

3  x  105 

2  x  106 

Total 

4  x  106 

4  X  106 

>  2  x  107 

Transition  layer  pressure  (dyne  cm-2) 

0.2 

0.07 

2 

Coronal  temperature  (K)  at  1.1  Rq 

1. 1-1.6  x  106 

106 

2.5  x  10° 

Coronal  energy  losses  (ergem-2  s-1) 

Conduction 

2  x  105 

6  x  104 

105  107 

Radiation 

105 

104 

5  x  106 

Solar  Wind 

<  5  x  104 

7  x  105 

<  105 

Total 

3  x  10s 

8  x  105 

107 

Solar  wind  mass  loss  (gcm~2s-1) 

<  2  x  10-11 

2  x  10-1° 

<  4  x  10-11 

43 


THERMONUCLEAR  FUSION 


26 


Natural  abundance  of  isotopes: 

hydrogen  nD/nH  =  1*5  *  10-4 

helium  nHe3 /nHe4  =  1*3  x  10-6 

lithium  nLi6/nLi7  =0-08 

Mass  ratios:  me/mD  =  2.72  x  10-4  =  1/3670 

(me/mD)1/2=  1.65  X  10-2  =  1/60.6 
mc/mT  =  1.82  x  10-4  =  1/5496 
(me/mT)1/2  =  1.35  x  10-2  =  1/74.1 

Absorbed  radiation  dose  is  measured  in  rads:  1  rad  =  102ergg-1.  The 
curie  (abbreviated  Ci)  is  a  measure  of  radioactivity:  1  curie  =  3.7  x 
1010  counts  sec-1. 

Fusion  reactions  (branching  ratios  are  correct  for  energies  near  the  cross 
section  peaks;  a  negative  yield  means  the  reaction  is  endothermic):27 


(la) 

D  +  D  50% 

♦T(l.OlMeV)  4-  p(3.02  MeV) 

(lb) 

50% 

+He3(0.82  MeV)  4-  n(2.45MeV) 

(2) 

D  +  T  - 

♦He4(3.5MeV)  4-  n(14.1MeV) 

(3) 

D  +  He3 - 

♦He4 (3.6  MeV)  4-  p(14.7MeV) 

(4) 

T  +  T  - 

♦He4  4-  2n  4-  11.3  MeV 

(5a) 

He3  +  T  51% 

♦  He4  4-  P  4-  n  4-  12.1  MeV 

(5b) 

43% 

♦He4(4.8  MeV)  4-  D(9.5  MeV) 

(5c) 

6% 

♦He5(2.4  MeV)  4-  p(U.9MeV) 

(6) 

P  +  Li6 

♦He4(1.7  MeV)  4-  He3(2.3MeV) 

(7a) 

P  +  Ll  20% 

♦2  He4  4-  17.3  MeV 

(7b) 

80% 

♦Be7  4-  n  —  1.6  MeV 

(8) 

D  4-  Li6  - 

♦2He4  4-  22.4  MeV 

(9) 

p  +  B11  - 

♦3  He4  4-  8.7  MeV 

(10) 

n  4-  Li6  - 

♦He4(2.1  MeV)  4-  T(2.7MeV) 

The  total  cross  section  in  ’  .ms  (1  barn  =  10-24  cm2)  as  a  function  of 
E,  the  energy  in  keV  of  the  incident  particle  [the  first  ion  on  the  left  side 
of  Eqs.  (l)-(5)j,  assuming  the  target  ion  at  rest,  can  be  fitted  by28 


(Tt(E)  = 


A5  4- 

(A4  —  A3E)2  4-  lj  A2 

E 

J^exp(Aj  E~1/7)  -  1  j 
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where  the  Duane  coefficients  Aj  for  the  principle  fusion  reactions  are  as  follows: 


1 

D-D 

(la) 

D-D 

(lb) 

D-T 

(2) 

D-He3 

(3) 

— 

T-He3 

(5a-c) 

i 

46.097 

47.88 

45.95 

89.27 

38.39 

■gw 

R 

372 

482 

50200 

25900 

448 

Qj 

4.36  x  10“4 

3.08  x  10-4 

1.368  x  10-2 

3.98  x  10-3 

1.02  x  10-3 

0 

sjj 

1.220 

1.177 

1.076 

1.297 

2.09 

0 

E 

0 

0 

409 

647 

0 

0 

Reaction  rates  av  (in  cm3  sec-1),  averaged  over  Maxwellian  distributions: 


Temperature 

(keV) 

D-D 

(la+  lb) 

1  <N 

Q  w 

D-He3 

(3) 

■9 

T-He3 

(5a-c) 

1.0 

1.5  x  10“22 

5.5  x  10-21 

10-26 

3.3  x  10-22 

10-28 

5.4  x  10“21 

2.6  x  10~19 

1.4  x  10-23 

7.1  x  10-21 

10-2S 

5.0 

1.8  X  10-19 

1.3  x  10“17 

6.7  X  10-21 

1.4  x  10-19 

2.1  x  10-22 

10.0 

1.2  x  10-18 

1.1  x  10“16 

2.3  x  10-19 

7.2  x  10-19 

1.2  x  10-2° 

20.0 

5.2  x  10-18 

4.2  x  10“16 

3.8  x  10“18 

2.5  x  10-18 

2.6  x  10-19 

50.0 

2.1  x  10“17 

8.7  x  10-16 

5.4  x  10-17 

8.7  x  10“18 

5.3  x  10-18 

100.0 

4.5  x  10-17 

8.5  x  10-16 

1.6  x  10-16 

1.9  x  10-17 

2.7  x  10-17 

200.0 

8.8  X  10-17 

6.3  x  10-16 

2.4  x  10-16 

4.2  x  10-17 

9.2  x  10-17 

500.0 

1.8  x  10-16 

3.7  x  10“16 

2.3  x  10-16 

8.4  x  10-17 

2.9  x  10-16 

1000.0 

2.2  X  10-16 

2.7  x  10"16 

1.8  x  10-16 

8.0  x  10-17 

5.2  x  10-16 

For  low  energies  (T  <  25  keV)  the  data  may  be  represented  by 

(vv)dd  -  2.33  x  10-14T-2/3  exp(-18.76T-1/3)  cm3  sec-1; 


(<tv)dt  —  3.68  x  10  12T  2^3  exp(  — 19.94 T  x^3)  cm3  sec  1 , 
where  T  is  measured  in  keV. 

The  power  density  released  in  the  form  of  charged  particles  is 

Pdd  —  3.3  X  10-  13ri£> 2 (Wv)dd  watt  cm ~ 3  (including  the  subsequent 

D-T  reaction); 

Pdt  =  5.6  x  10-13n£>n'r(ol7)£>r  watt  cm-3; 

P 0He^  —  2.9  X  10  Watt  CHI 
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RELATIVIS  fIC  ELECTRON  BEAMS 


Here  7  =  (1  —  /32)-1/2  is  the  relativistic  scaling  factor;  quantities 
in  analytic  formulas  are  expressed  in  SI  or  cgs  units,  as  indicated;  in 
numerical  formulas,  I  is  in  amperes  (A),  B  is  in  gauss  (G),  electron 
linear  density  N  is  in  cm-1,  and  temperature,  voltage  and  energy  are  in 
MeV;  /3*  =  vzf  c;  fc  is  Boltzmann’s  constant. 

Relativistic  electron  gyroradius: 

(72  -  1)1/2  (cgs)  =  1.70  x  103(72  -  l)1/2£-1  cm. 
Relativistic  electron  energy: 


W  =  me2 7  =  0.5 1 17  MeV. 


Bennett  pinch  condition: 

/ 2  =  2 Nk(Tt  +  T,)c2  (cgs)  =  3.20  x  10_4AT(TC  +  T{)  A2. 
Alfven-Lawson  limit: 

IA  =  (mc3/e)B, 7  (cgs)  =  (47TTOc//i0e)/?*7  (SI)  =  1.70  x  10*/3zjA. 
The  ratio  of  net  current  to  IA  is 


J_  _  v_ 
fA  7  ’ 

Here  v  —  Nrc  is  the  Budker  number,  where  re  =  e2/mc2  =  2.82  x 
10  cm  is  the  classical  electron  radius.  Beam  electron  number  density 
is 

nb  =  2.08  x  10 cm-3, 

where  J  is  the  cuirent  density  in  A  cm-2.  For  a  uniform  beam  of  radius 
a  (in  cm), 

rib  =  G.63  x  10‘  Ia~2 1  cm-3, 


an  d 


2  Tf.  v 

a  7 
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Child's  law:  (nun-relativistic)  space-charge-limited  current,  density  be¬ 
tween  [>arallel  jdates  with  voltage  drop  V  and  separation  d  (in  cm)  is 

J  =  2.34  x  10 3V3/2d~2  A  cm-2. 


The  saturated  parapotential  current  (magnetically  self-limited  flow  along 
equipotentials  in  pinched  diodes  and  transmission  lines)  is29 


Ip  —  8.5  x  103G7ln  7  -f  (j2  —  l)1^2  A, 


where  G  is  a  geometrical  factor  depending  on  the  diode  structure: 

w  for  parallel  plane  cathode  and  anode 

2-rcd  of  width  w,  separation  d ; 

/  B,\-' 

G 


G  = 


■  (•■  £) 


for  cylinders  of  radii  R i  (inner)  and  Rz  (outer); 


For  P  — *  0  (7 


for  conical  cathode  of  radius  Rc ,  maximum 
separation  do  (at  r  =  Rc)  from  plane  anode. 
1),  both  I  a.  and  Ip  vanish. 


The  condition  for  a  longitudinal  magnetic  field  Bz  to  suppress  filamen- 
tation  in  a  beam  of  current  density  J  (in  A  cm-2)  is 

Bt  >  47/M7J)1/2G. 


Voltage  registered  by  Rogowski  coil  of  minor  cross-sectional  area  A,  tl 
turns,  major  radius  a,  inductance  L ,  external  resistance  R  and  capaci¬ 
tance  C  (all  in  SI): 


externally  integrated 
self-integrating 


V  =  (l/7?C)(nA/*0//2ira); 

V  —  (R/ L)(nAfio  I/2ira)  =  Rl/n. 


X-ray  production,  target  with  average  atomic  number  Z  ( V  <  5  MeV): 
77  =  x-ray  power/beam  power  =  7  x  10-4ZV. 

X-ray  dose  at  I  meter  generated  by  an  e-beam  depositing  total  charge  Q 
coulombs  while  V  >  0.84Vmax  >n  material  with  charge  state  Z : 


D  =  150 Vml  *QZ1/2  rads. 
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BEAM  INSTABILITIES30 


Name 

Conditions 

Saturation  Mechanism 

Electron- 

electron 

Vd  >  Vei ,  j  1,2 

Electron  trapping  until 

Vej  ~  Vd 

Buneman 

Vd  >  ( M  j m)1/3  Vi , 
Vd  >  K 

Electron  trapping  until 

K  ~  Vd 

Beam-plasma 

Vb  >  {np/nh)x'zVb 

Trapping  of  beam  electrons 

Weak  beam- 
plasma 

Vb  <  (  Tlp  / rif,)1/3  Vb 

Quasilinear  or  nonlinear 
(mode  coupling) 

Beam-plasma 

(hot-electron) 

ve  >  Vb  >  Vb 

Quasilinear  or  nonlinear 

Ion  acoustic 

Te  >  T, ,  Vd  >  C. 

Quasilinear,  ion  tail  form¬ 
ation,  nonlinear  scattering, 
or  resonance  broadening. 

Anisotropic 

temperature 

(hydro) 

Tex  >  2Tej{ 

Isotropization 

Ion  cyclotron 

Vd  >  20 Vi  (for 

Te  «  Ti) 

Ion  heating 

Beam-cyclotron 

(hydro) 

vd  >  Ca 

Resonance  broadening 

Modified  two- 
stream  (hydro) 

Vd  <  (1  +  (3)l/2VA , 
Vd  >  c. 

Trapping 

Ion-ion  (equal 
beams) 

U  <  2(1  +  0)l/2VA 

Ion  trapping 

Ion-ion  (equal 
beams) 

U  <  2Ca 

Ion  trapping 

For  nomenclature,  see  p.  50. 
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Parameters  of  Most  Unstable  Mode 


Name 

Wave 

Group 

Growth  Rate 

Frequency 

Number 

Velocity 

1 

UJ, 

Electron- 

-UJe 

0 

0.9  — 

0 

electron 

2 

Vd 

Buneman 

(  m  \  1/3 

u>e 

2 

0.7  (  —  )  u>e 

0.4  (  —  J  u>e 

-Vd 

\M  J 

Vd 

3 

Beam-plasma 

0.7  (  —  )  u>e 

uc  - 

2 

-Vb 

\nP  J 

nJn<-Y'3 

Vb 

3 

0.4  I  -  1  u>e 

\nP  J 

Weak  beam- 

UJe 

3Ve2 

plasma 

2 np  \VbJ 

vb 

v6 

Beam-plasma 

(hot-electron) 

\nj  Vb 

vb 

Ve 

v6 

{  m\1/2 

Ion  acoustic 

U) 

Ul, 

ca 

Anisotropic 

u>e  cos  6  ~  fte 

r7l 

temperature 

(hydro) 

- 

1  - 
3Vi 

Ion  cyclotron 

0.1ft/ 

1.2  ft/ 

r"1 

1 

Beam-cyclotron 

0.7fte 

7lfte 

0.7A-1 

£  Vd\ 

(hydro) 

—  ft H 

<  c8 

Modified  two- 

0.9  ftH 

ft// 

1.7 

1 

-vd 

stream 

2 

Vd 

2 

(hydro) 

ft// 

1.2 

U 

Ion-ion  (equal 

0.4ftw 

0 

0 

beams) 

Ion-ion  (equal 

0.4a;, 

0 

u >, 

1.2  — 

0 

beams) 

U 

For  nomenclature,  see  p.  50. 
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In  the  preceding  tables,  subscripts  e,  d,  b,  p  stand  for  “electron,” 
“ion,”  “drift,”  “beam,”  and  “plasma,”  respectively.  Thermal  velocities 
arc  denoted  by  a  bar.  In  addition,  the  following  are  used: 


m 

electron  mass 

re,n 

gyroradius 

M 

ion  mass 

(1 

plasma/magnetic  energy 

V 

velocity 

density  ratio 

T 

temperature 

Va 

Alfven  speed 

number  density 

Qe,Qi 

gyrofrequency 

n 

harmonic  number 

Qh 

hybrid  gyrofrequency, 

Ca  =  ( TJM )1/2 

ion  sound  speed 

Qh2  =  c  G  i 

,  t J, 

plasma  frequency 

U 

relative  drift  velocity  of 

Debye  length 

two  ion  species 

LASERS 


System  Parameters 

Efficiencies  and  power  levels  are  approximately  state-of-the-art  (1990). 31 


Type 

Wavelength 

Efficiency 

Power  levels  available  (W) 

(/im) 

Pulsed 

cw 

co2 

10.6 

0.01-0.02 

(pulsed) 

>  2  X  1013 

>  105 

CO 

5 

0.4 

>  109 

>  100 

Holmium 

2.06 

0.031-0. It 

>  107 

30 

Iodine 

1.315 

0.003 

>  1012 

- 

Nd-glass, 

YAG 

1.06 

0.001-0.061 
>  o.it 

~  1014  (ten- 
beam  system) 

1-103 

♦Color  ceuter 

1-4 

10-3 

>  106 

1 

♦Vibronic  (Ti 
Sapphire) 

0.7  0.9 

>  0.1  x  i)p 

10G 

1-5 

Ruby 

0.6943 

<  10-3 

O 

O 

f-H 

1 

He-Ne 

0.6328 

10~4 

1  50xl0~3 

♦Argon  ion 

0.45-0.60 

10~3 

5  x  104 

1  20 

*OPO 

0. 4-9.0 

>  0.1  X  T]p 

106 

15 

n2 

0.3371 

0.001  -0.05 

105-1Gc 

♦Dye 

0.3  1.1 

10~s 

>  10c 

140 

Kr-F 

0.26 

0.08 

>  10° 

500 

Xenon 

0.175 

0.02 

>  10s 

♦Tunable  sources  jlamp-driven  Jdiode-drivcn 

YAG  stands  for  Yttrium  Aluminum  Garnet  and  OPO  for  Optical  Para¬ 
metric  Oscillator;  rjp  is  pump  laser  efficiency. 
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Formulas 

An  e-m  wave  with  k  j|  B  has  an  index  of  refraction  given  by 
7*±  =  [1  -  ^  a»ce)]1/2, 

where  ±  refers  to  the  helicity.  The  rate  of  change  of  polarization  angle 
9  as  a  function  of  displacement  a  (Faraday  rotation)  is  given  by 

dO/ds  =  ( fc/2)(n _  -  n+)  =  2.36  x  10 *NBf~2  cm-1, 

where  N  is  the  electron  number  density,  D  is  the  field  strength,  and  /  is 
the  wave  frequency,  all  in  cgs. 

The  quiver  velocity  of  an  electron  in  an  e-m  field  of  angular  fre¬ 
quency  u>  is 


vo  =  e  FInax  /  rnu>  =  25.6/1/,2Ao  cm  sec  1 

in  terms  of  the  laser  flux  I  =  cFmax/ 8ir,  with  I  in  watt/cm2,  laser 
wavelength  A0  in  /im.  The  ratio  of  quiver  energy  to  thermal  energy  is 

Wqu/Wth  =  mev02/2kT  =  1.81  x  10_13A02//T, 

where  T  is  given  in  eV.  For  example,  if  I  —  1015  Wcm-2,  A0  =  1  /im, 
T  =  2  keV,  then  Wqu/VFlh  %  0.1. 

Pondermotive  force: 

T  =  NV(E2)/8*Nc, 

where 

Nc  =  1.1  x  1021  A0  -2cin-3. 

For  uniform  illumination  of  a  lens  with  /-number  F,  the  diameter 
d  at  focus  (85%  of  the  energy)  and  the  depth  of  focus  l  (distance  to  first 
zero  in  intensity)  are  given  by 

d  zz  2A4F\9/6ol  and  l  %  -F.2F2 \9  j 9 ol  ■ 

Here  9  is  the  beam  divergence  containing  85%  of  energy  and  9dl  is  the 
diffraction -limited  divergence: 

9  dl  —  2.44A/6, 

where  b  is  the  aperture.  These  formulas  are  modified  for  nonuniform 
(such  .as  Gaussian)  illumination  of  the  lens  or  for  pathological  laser  pro¬ 
files. 
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ATOMIC  PHYSICS  AND  RADIATION 


Energies  and  temperatures  are  in  eV;  all  other  units  are  cgs  except 
where  noted.  Z  is  the  charge  state  (Z  =  0  refers  to  a  neutral  atom); 
the  subscript  e  labels  electrons.  N  refers  to  number  density,  n  to  princi¬ 
pal  quantum  number.  Asterisk  superscripts  on  level  population  densities 
denote  local  thermodynamic  equilibrium  (LTE)  values.  Thus  Nn*  is  the 
LTE  number  density  of  atoms  (or  ions)  in  level  n. 

Characteristic  atomic  collision  cross  section: 


7rao2  =  8.80  X  10  17  cm2. 


Binding  energy  of  outer  electron  in  level  labelled  by  quantum  numbers 
n,  /: 


EL(nJ)  =  ~ 


z2eh 

oo 

(n  -  A,)2’ 


where  E ^  =  13.6  eV  is  the  hydrogen  ionization  energy  and  Aj  =  0.75/  5, 
1  >  5,  is  the  quantum  defect. 

Excitation  and  Decay 

Cross  section  (Bethe  approximation)  for  electron  excitation  by  dipole 
allowed  transition  m  —*  n  (Refs.  32,  33): 

—  13  f  n  m  Q  (  ^  i  )  2 

(3)  o-nn  =  2.36  x  10  - — -  cm  , 

eA  Enru 

where  /nm  >s  the  oscillator  strength,  g(n,m)  is  the  Gaunt  factor,  e  is  the 
incident  electron  energy,  and  A En  m  —  En  E  rlx  . 

Electron  excitation  rate  averaged  over  Maxwellian  velocity  distribution, 
Xmn  =  (<Tm„i>)  (Refs.  34,  35): 


-5  fn  m  <tf(n,  m))Nr  (  A En 

(4)  Xmn=  1.6x10  - - - —jz  exp  ( - — 

A£„  rl 7  2  V 


m  \  —  1 

—  sec 


where  ( g(n ,  m))  denotes  the  thermal  averaged  Gaunt  factor  (generally  ~  1 
for  atoms.  ~  0.2  for  ions). 
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Hate  for  electron  eollisional  deexcitation: 


(5)  Ynm  =  (Nru  */Nn*)X . . 

Here  N m* / N n*  —  (3m/3n)cxp(A£nm/Tc)  is  the  Boltzmann  relation  foi 
level  population  densities,  where  gn  is  the  statistical  weight  of  level  n. 

Rate  for  spontaneous  decay  n  — *  m  (Einstein  A  coefficient)34 

(6)  i4nm  =  4.3  x  107(3n/5m)/r)m(A£nm)2sec"1. 


Intensity  emitted  per  unit  volume  from  the  transition  n  — ♦  m  in  an  optic¬ 
ally  thin  plasma: 

(7)  Inrn  =  16  x  1 0  ~  1 9  A  n  m  Nn  A  En  m  watt/cm3. 


Condition  for  steady  state  in  a  corona  model: 


(8) 


N0Ne(a0nv)  =  Nn  Ano, 


where  the  ground  state  is  labelled  by  a  zero  subscript. 

Hence  for  a  transition  n  — ►  m  in  ions,  where  {<7(71,  0))  0.2, 


(9) 


=  5.1  x  10 


—  25  fnm  9 0 


NeN0  / 

a/2  l 


A  Er 


A  E 


exp 


n  0 


AJS'uo  \  watt 
Tc  /  cm3 


Ionization  and  Recombination 

In  a  general  time-dependent  situation  the  number  density  of  the 
charge  state  Z  satisfies 


(10) 


dN(Z) 

dt 


= 


-  S(Z)N{  Z)  -  a(Z) N(Z) 


+  S(Z  -  1)N(Z  —  1)  +  a(Z  +  l)N(Z  +  1) 


Here  S(Z)  is  the  ionization  rate.  The  recombination  rate  c*(Z)  has  the 
form  c*{Z)  =  ar{Z)  4-  Nfaj(Z),  where  ar  and  0:3  are  the  radiative  and 
three-body  recombination  rates,  respectively. 
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Classical  ionization  cross-section36  for  any  atomic  shell  j 


(11)  <t,  =  G  x  10  14  bj  g  j(x)  /  U  j2  cm2 . 

Here  b:  is  the  number  of  shell  electrons;  Uj  is  the  binding  energy  of  the 
ejected  electron;  x  —  e/U},  where  e  is  the  incident  electron  energy;  and  g 
is  a  universal  function  with  a  minimum  value  gmin  «  0.2  at  i  «  4. 

Ionization  from  ion  ground  state,  averaged  over  Maxwellian  electron  dis¬ 
tribution,  for  0.02  <  Te/E ^  <  100  (Ref.  35): 


(12) 


S(Z )  =  10 


-5 


(TJEi) 


1/2 


(££)3/2(6.0  +  TJEZ) 


exp 


where  E ^  is  the  ionization  energy. 

Electron-ion  radiative  recombination  rate  (e  -f  N(Z)  — ►  N(Z  —  1)  -f  hu) 
for  TJZ2  <  400eV  (Ref.  37): 


(13)  ar(Z)  =  5.2  x  10 


-14 


(#) 


Z  \  1/2 


0.43+  ill.  (E^/T,) 


+0.469(B^,/Te)_1/3 


3  / 

cm  /sec. 


For  1  eV  <  Te  / Z2  <  15  eV,  this  becomes  approximately35 

(14)  ar(Z)  =  2.7  x  10~ 13 Z2 Tc  ~  1/^2  cm3/scc. 

Collisional  (three-body)  recombination  rate  for  singly  ionized  plasma:38 

(15)  «3  =  8.75  x  10~27  Te  ~4  5  cm6 /sec. 


Photoionization  cross  section  for  ions  in  level  n,  l  (short- wavelength  limit): 

(16)  ap\t(n,l)  =  1.64  x  10_  10  Z5/n3  K  ‘  cm2 . 

where  K  is  the  wavenumber  in  Rydbergs  ( 1  Rydberg  =  1 .0974  X  10°  cm -  1  ). 
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—  3 


Ionization  Equilibrium  Models 


Saha  equilibrium: 


39 


(17) 


Ne  Ni*(Z) 
Nn*(Z  ~  1) 


g?Te  3/2  /  Ez(n,l)\ 

=  6.0  x  102  - ; —  exp  ( - 2^ - J  cm 


Z-  1 


where  gZ  is  the  statistical  weight  for  level  n  of  charge  state  Z  and  £^,(n,  l) 
is  the  ionization  energy  of  the  neutral  atom  initially  in  level  (n.,1),  given 
by  Eq.  (2). 

In  a  steady  state  at  high  electron  density, 


(18) 


NCN*(Z)  S(Z  -  1) 


N*(Z  -  1) 


«3 


a  function  only  of  T. 

Conditions  for  LTE:39 

(a)  Collisional  and  radiative  excitation  rates  for  a  level  n  must  satisfy 


(19) 


Yn  m  ~  1 0  .A  n  m  ■ 


(b)  Electron  density  must  satisfy 


(20) 


Ne  >  7  x  1018  Z7  n~  17 /7  (T/  E^)1/2  cm'3 


Steady  state  condition  in  corona  model: 


(21) 


N(Z  -  1) 
N(Z)~ 


a. 


S(Z  -  1) 


Corona  model  is  applicable  if40 


(22) 


10wt,  <  Ne  <  I0lb  Te  ‘ ~  cm  ~  J  , 


where  t  /  is  the  ionization  time. 
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Ftadiat  ion 


N.  D.  Energies  and  temperatures  are  in  eV;  all  other  quantities  are 
in  cgs  units  except  where  noted.  Z  is  the  charge  state  (Z  =  0  refers  to  a 
neutral  atom);  the  subscript  e  labels  electrons.  N  is  number  density. 

Average  radiative  decay  rate  of  a  state  with  principal  quantum  number  n 
is 

(23)  An  =  ^  Anm  =  1.6  X  10 1QZ4n~9/2  sec. 

m  <  n 


Natural  linewidth  (A E  in  eV): 

(24)  A£  At  =  h  =  4.14  x  10-15  eVsec, 

where  At  is  the  lifetime  of  the  line. 

Doppler  width: 

(25)  AA/A  =  7.7  x  10 “5(T//i)1/2, 

where  n  is  the  mass  of  the  emitting  atom  or  ion  scaled  by  the  proton  mass. 
Optical  depth  for  a  Doppler-broadened  line:39 

(26)  r  =  3.52  x  10_13/nm  A( Me2/ fcT) 1/2  NL  =  5.4  x  10-9 X(n/T)l/2  NL, 

where  /„m  is  the  absorption  oscillator  strength,  A  is  the  wavelength,  and 
L  is  the  physical  depth  of  the  plasma;  M,  N ,  and  T  are  the  mass,  number 
density,  and  temperature  of  the  absorber;  ft  is  M  divided  by  the  proton 
mass.  Optically  thin  means  r  <  1. 

Resonance  absorption  cross  section  at  center  of  line: 

(27)  cr\  =  \c  —  5.6  X  10- 13  A2/ A  A  cm2 . 

Wien  displacement  law  (wavelength  of  maximum  black-body  emission): 

(28)  Amax  =  2.50  x  10~5r_l  cm. 

Radiation  from  the  surface  of  a  black  body  at  temperature  T: 

(29)  W  =  1.03  x  10ST4  watt/cm2. 
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BremsstraMung  from  hydrogcn-like  plasma:26 


(30) 


PB r  =  1.69  x  10~32A^cTe1/2  ^2  [z7N{Z)^  watt/ 


cm 


where  the  sum  is  over  all  ionization  states  Z. 
Bremsstrahlung  optical  depth:41 


(31)  r  =  5.0  x  10"38JVC  N,Z2gLT~7/2, 

where  ~g  «  1.2  is  an  average  Gaunt  factor  and  L  is  the  physical  path  length. 

Inverse  bremsstrahlung  absorption  coefficient42  for  radiation  of  angular 
frequency  u>: 


(32)  x  =  3.1  x  10-7Zne2  In  AT'3/2w"2(l  -  u>2p/u2)1/2  cm-1; 

here  A  is  the  electron  thermal  velocity  divided  by  V,  where  V  is  the  larger 
of  u>  and  u>p  multiplied  by  the  larger  of  Ze2  /kT  and  h/(mkT)1^2 

Recombination  (free-bound)  radiation: 


(33) 


Pr  =  1.69  X  10 


-32 


Z2N(Z) ^ 


E 


z- 1 


watt/cm* 


Cyclotron  radiation26  in  magnetic  field  B: 

(34)  Pc  =  6  21  x  10~2S  B2  NeTe  watt/cm3. 

For  NekTe  =  NikTt  =  D2  /  I6n  ((3  =  1,  isothermal  plasma),26 

(35)  Pc  =  5.00  x  10-387V2T2  watt/cm3. 


Cyclotron  radiation  energy  loss  e-folding  time  for  a  single  electron:41 

9.0  x  108P-2 


(36) 


tc 


2.5  +  7 


sec. 


where  7  is  the  kinetic  plus  rest  energy  divided  by  the  rest  energy  me2. 
Number  of  cyclotron  harmonics41  trapped  in  a  medium  of  finite  depth  L: 

(37)  mtr  =  (57/?PL)1/6, 


where  ,3  =  8tt  N  kT /  D2  . 

Line  radiation  is  given  by  summing  Eq.  (9)  over  all  species  in  the  plasma. 
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ATOMIC  SPECTROSCOPY 


Spectroscopic  notation  combines  observational  and  theoretical  el¬ 
ements.  Observationally,  spectral  lines  are  grouped  in  series  with  line 
spacings  which  decrease  toward  the  series  limit.  Every  line  can  be  related 
theoretically  to  a  transition  between  two  atomic  states,  each  identified  by 
its  quantum  numbers. 

Ionization  levels  are  indicated  by  roman  numerals.  Thus  C  I  is  un¬ 
ionized  carbon,  CII  is  singly  ionized,  etc.  The  state  of  a  one-electron 
atom  (hydrogen)  or  ion  (Hell,  Li  III,  etc.)  is  specified  by  identifying  the 
principal  quantum  number  n  =  1,2,...,  the  orbital  angular  momentum 
/  =  0,1,...,  n  —  1,  and  the  spin  angular  momentum  s  =  i|.  The  total 
angular  momentum  j  is  the  magnitude  of  the  vector  sum  of  1  and  s,  j  = 
l  ±  5  (j  >  j).  The  letters  s,  p,  d,  f,  g,  h,  i,  k,  1,  ...,  respectively,  are 
associated  with  angular  momenta  l  =  0,  1,  2,  3,  4,  5,  6,  7,  8,  ...  .  The 
atomic  states  of  hydrogen  and  hydrogenic  ions  are  degenerate:  neglecting 
fine  structure,  their  energies  depend  only  on  n  according  to 

_  R00hcZ2n~2  _  RyZ2 
n  1  -f  m/M  n2 

where  h  is  Planck’s  constant,  c  is  the  velocity  of  light,  m  is  the  electron 
mass,  M  and  Z  are  the  mass  and  charge  state  of  the  nucleus,  and 

Roc  =  109,737  cm-1 

is  the  Rydberg  constant.  If  En  is  divided  by  he,  the  result  is  in  wavenum¬ 
ber  units.  The  energy  associated  with  a  transition  m  — ►  n  is  given  by 

A  £mn  =  Ry(  1  /  tji2  -  1/n2), 

with  m  <  n  (m  >  n)  for  absorption  (emission)  lines. 

For  hydrogen  and  hydrogenic  ions  the  series  of  lines  belonging  to  the 
transitions  rn  — *  n  have  conventional  names: 


Transition 

1  — ♦  n 

2  — »  n 

3  — *  71 

4  — *  Tl 

5  — *  71 

*0 

t- 

I 

0 

Name 

Lyman 

Balmer 

Paschen 

Brackett 

Pfund 

Humphreys 

Successive  lines  in  any  series  are  denoted  nr,  /?,  7,  etc.  Thus  the  transi¬ 
tion  1  — > >  3  gives  rise  to  the  Lyman-/?  line.  Relativistic  effects,  quantum 
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electrodynamic  effects  (e.g.,  the  Lainb  shift),  and  interactions  between 
the  unclear  magnetic  moment  and  the  magnetic  field  due  to  the  electron 
produce  small  shifts  and  splittings  (  <  10-2cm-1);  the  last  of  these  are 
called  “hyperfine  structure.” 

In  many-electron  atoms  the  electrons  are  grouped  in  closed  and  open 
shells,  with  spectroscopic  properties  determined  mainly  by  the  outer  shell. 
Shell  energies  depend  primarily  on  n;  the  shells  corresponding  to  n  =  1,  2, 
3, .  .  .  are  called  K ,  L,  M ,  etc.  A  shell  is  made  up  of  subshells  of  different 
angular  momenta,  each  labeled  according  to  the  values  of  n,  l ,  and  the 
number  of  electrons  it  contains  out  of  the  maximum  possible  number, 
2(2/  +  1).  For  example,  2p5  indicates  that  there  are  5  electrons  in  the 
subshell  corresponding  to  /  =  1  (denoted  by  p)  and  n  —  2. 

In  the  lighter  elements  the  electrons  fill  up  subshells  within  each  shell 
in  the  order  s,  p,  d,  etc.,  and  no  shell  acquires  electrons  until  the  lower 
shells  are  full.  In  the  heavier  elements  this  rule  does  not  always  hold.  But 
if  a  particular  subshell  is  filled  in  a  noble  gas,  then  the  same  subshell  is 
filled  in  the  atoms  of  all  elements  that  come  later  in  the  periodic  table. 
The  ground  state  configurations  of  the  noble  gases  are  as  follows: 

He  Is2 

Ne  ls22s22p6 

Ar  ls22s22p63s23p6 

Kr  ls22s22p63s23p63d104s24p6 

Xe  ls22s22p63s23p63d104s24p64d105s25p6 

Rn  ls22s22p63s23p63d104s24p64d104f  145s25p65d106s26p6 

Alkali  metals  (Li,  Na,  K,  etc.)  resemble  hydrogen;  their  transitions  are 
described  by  giving  n  and  /  in  the  initial  and  final  states  for  the  single 
outer  (valence)  electron. 

For  general  transitions  in  most  atoms  the  atomic  states  are  specified 
in  terms  of  the  parity  (  —  1)E,»  and  the  magnitudes  of  the  orbital  angular 
momentum  L  =  El,-,  the  spin  S  =  Es,,  and  the  total  angular  momentum 
J  =  L  +  S,  where  all  sums  are  carried  out  over  the  unfilled  subshells  (the 
filled  ones  sum  to  zero).  If  a  magnetic  field  is  present  the  projections 
Ms,  and  M  of  L,  S,  and  J  along  the  field  are  also  needed.  The  quantum 
numbers  satisfy  |M/,|  <  L  <  vl\  |Ms|  <  S  <  v/ 2;  and  \M\  <  J  <  L  +  S , 
where  v  is  the  number  of  electrons  in  the  unfilled  subshell.  Upper  case 
letters  S,  P,  D,  etc.,  stand  for  L  =  0,1,  2,  etc.,  in  analogy  with  the 
notation  for  a  single  electron.  For  example,  the  ground  state  of  Cl  is 
described  by  3p5  2P3/2-  The  first  part  indicates  that  there  are  5  electrons 
in  the  subshell  corresponding  to  n  =  3  and  /  =  1 .  (The  closed  inner  sub¬ 
shells  ls22s22p63s2,  identical  with  the  configuration  of  Mg,  are  usually 
omitted.)  The  symbol  lP’  indicates  that  the  angular  momenta  of  the  outer 


59 


electrons  combine  to  give  L  =  1.  The  prefix  ‘2’  represents  the  value  of  the 
multiplicity  2 S  +  1  (the  number  of  states  with  nearly  the  same  energy), 
whi  ch  is  equivalent  to  specifying  S  —  ^ .  The  subscript  3/2  is  the  value  of 
,7.  The  superscript  ‘o’  indicates  that  the  state  has  odd  parity;  it  would  be 
omitted  if  the  state  were  even. 

The  notation  for  excited  states  is  similar.  For  example,  helium  has  a 
state  ls2s  3S;  which  lies  19.72eV  (159,856cm-1)  above  the  ground  state 
Is2  1  So-  But  the  two  “terms”  do  not  “combine”  (transitions  between  them 
do  not  occur)  because  this  would  violate,  e.g.,  the  quantum-mechanical 
selection  rule  that  the  parity  must  change  from  odd  to  even  or  from  even 
to  odd.  For  electric  dipole  transitions  (the  only  ones  possible  in  the  long- 
wavelength  limit),  other  selection  rules  are  that  the  value  of  7  of  only  one 
electron  can  change,  and  only  by  A 7  =  ±1;  AS  =  0;  A L  =  ±1  or  0;  and 
A  J  —  ±1  or  0  (but  L  =  0  does  not  combine  with  L  =  0  and  J  =  0  does  not 
combine  with  J  =  0).  Transitions  are  possible  between  the  helium  ground 
state  (which  has  S  =  0,  L  =  0,  J  =  0,  and  even  parity)  and,  e.g.,  the 
state  ls2p  1  P°  (with  5  =  0,  L  =  1,  J  =  1,  odd  parity,  excitation  energy 
21.22  eV).  These  rules  hold  accurately  only  for  light  atoms  in  the  absence 
of  strong  electric  or  magnetic  fields.  Transitions  that  obey  the  selection 
rules  are  called  “allowed”;  those  that  do  not  are  called  “forbidden.” 

The  amount  of  information  needed  to  adequately  characterize  ?  state 
increases  with  the  number  of  electrons;  this  is  reflected  in  the  notation. 
Thus43  Oil  has  an  allowed  transition  between  the  states  2p23p'  2F2/,2 

and  2p2(3D)3d'  2F7/2  (and  between  the  states  obtained  by  changing  J 
from  7/2  to  5/2  in  either  or  both  terms).  Here  both  states  have  two 
electrons  with  n  =  2  and  7  =  1;  the  closed  subshells  ls22s2  are  not  shown. 
The  outer  (n.  =  3)  electron  has  7  =  1  in  the  first  state  and  7  =  2  in  the 
second.  The  prime  indicates  that  if  the  outermost  electron  were  removed 
by  ionization,  the  resulting  ion  would  not  be  in  its  lowest  energy  state. 
The  expression  (1D)  give  the  multiplicity  and  total  angular  momentum  of 
the  “parent”  term,  i.e.,  the  subshell  immediately  below  the  valence  sub¬ 
shell;  this  is  understood  to  be  the  same  in  both  states.  (Grandparents, 
etc.,  sometimes  have  to  be  specified  in  heavier  atoms  and  ions.)  Another 
example43  is  the  allowed  transition  from  2p2(3P)3p  2P°^2  (or  2P^2)  to 

2p2(JD)3d'  2 S i / 2 1  in  which  there  is  a  “spin  flip”  (from  antiparallel  to 
parallel)  in  the  n  —  2,7  =  1  subshell,  as  well  as  changes  from  one  state  to 
the  other  in  the  value  of  7  for  the  valence  electron  and  in  L. 

The  description  of  fine  structure,  Stark  and  Zeeman  effects,  spectra 
of  highly  ionized  or  heavy  atoms,  etc.,  is  more  complicated.  The  most 
important  difference  between  optical  and  X-ray  spectra  is  that  the  latter 
involve  energy  changes  of  the  inner  electrons  rather  than  the  outer  ones; 
often  several  electrons  are  involved. 
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